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ABSTRACT 

We have studied a sample of 809 Mg II absorption systems with 1.0^ z a b s ^1.86 in the 
spectra of Sloan Digital Sky Survey QSOs, with the aim of understanding the nature and 
abundance of the dust and the chemical abundances in the intervening absorbers. Normalized, 
composite spectra were derived, for abundance measurements, for the full sample and several 
sub-samples, chosen on the basis of the line strengths and other absorber and QSO properties. 
Average extinction curves were obtained for the sub-samples by comparing their geometric 
mean spectra with those of matching samples of QSOs without absorbers in their spectra. 
There is clear evidence for the presence of dust in the intervening absorbers. The 2175 A 
feature is not present in the extinction curves, for any of the sub-samples. The extinction 
curves are similar to the SMC extinction curve with a rising UV extinction below 2200 A. The 
absorber rest frame colour excess, E(B — V), derived from the extinction curves, depends on 
the absorber properties and ranges from <0.001 to 0.085 for various sub-samples. The column 
densities of Mg II, Al II, Si II, Ca II, Ti II, Cr II, Mn II, Fe II, Co II, Ni II, and Zn II do not show 
such a correspondingly large variation. The overall depletions in the high E(B — V) samples 
are consistent with those found for individual damped Lyman alpha systems, the depletion 
pattern being similar to halo clouds in the Galaxy. Assuming an SMC gas-to-dust ratio we 
find a trend of increasing abundance with decreasing extinction; systems with Nh i ~ 10 20 
cm -2 show solar abundance of Zn. The large velocity spread of strong Mg II systems seems 
to be mimicked by weak lines of other elements. The ionization of the absorbers, in general 
appears to be low: the ratio of column densities of Al III to Al II is always less than 1/2. QSOs 
with absorbers are, in general, at least three times as likely to have highly reddened spectra as 
compared to QSOs without any absorption systems in their spectra. 
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1 INTRODUCTION 

1.1 The nature of QSO absorbers 

Absorption lines in spectra of quasi-stellar objects have been 
known since shortly after the discovery of QSOs (Sandage 1965; 
Schmidt 1966; Arp, Bolton & Kinman 1967; Burbidge, Lynds & 
Burbidge 1966; Stockton & Lynds 1966). While many of the first 
detections were close to the redshifts of the QSOs, what are now 
termed associated systems, Bahcall, Peterson & Schmidt (1966) 
found one of the earliest systems which could be termed interven- 
ing and Bahcall (1968) established the reality of such systems. Bah- 
call and Spitzer (1969) presented plausibility arguments that fore- 
ground galaxies with extended halos could produce some of the 
intervening systems. If such were the origin of the QSO absorption 
line systems (QSOALSs), they would provide an excellent way to 
understand various aspects of structure formation and evolution in 
the universe, in particular, the buildup of the elements over cos- 
mic time and evolution of element flows into and out of galaxies 
as they form (Hermann et al. 2001; Adelberger et al. 2003, 2005; 
Tumlinson & Fang 2005). 

There are, indeed, specific examples of QSO absorption lines 
coming from specific galaxies (Cohen et al. 1987; Yanny, York & 
Gallagher 1989; Bergeron 1988; Chen, Kennicut & Rauch 2005). 
On the other hand, there are numerous cases in which a QSO sight- 
line passes near a visible galaxy without producing absorption lines 
(Bechtold & Ellingson 1992). Some of the galaxies are almost cer- 
tainly dwarf galaxies (York et al. 1986; Rao et al. 2003; Haehnelt, 
Stenimetz & Rauch 1998). Others may have an outer morphology 
that resembles the Milky Way (MW) distribution of 21 cm clouds 
(York 1982; Wakker & van Woerden 1997). Such a patchy con- 
figuration would explain why absorbers are sometimes found near 
galaxies and why some QSO/near galaxy coincidences do not pro- 
duce absorbers. In some cases extended, turbulent regions around 
star-forming galaxies (Adelberger et al. 2005) may show up as 
QSOALSs. 

Despite numerous observations of absorption line systems 
in spectra of QSOs, there may be some selection effects in the 
samples. Why, for instance are QSOALSs with fully or near so- 
lar abundances seldom seen (Prochaska et al. 2003a; Khare et al. 
2004; Kulkarni et al. 2005a)? Only a few systems with metallic- 
ity above solar are definitely known: SDSSJ1323-0021, z a6s =0.716 
(Khare et al. 2004), Q0058+019, z a6s =0.612 (Pettini et al. 2000). 
There may be some selection effects that keep such absorbers out 
of the observed sample. For example, dwarf galaxies may consti- 
tute a large part of the sample; these usually have sub-solar abun- 
dances. Also, absorbers with solar abundances may have more dust, 
making any background objects fainter and keeping those QSOs 
out of magnitude-limited samples (e.g. Fall & Pei 1993; Boisse 
et al. 1998; Vladilo & Peroux 2005). Ostriker & Heisler (1984) 
argued that foreground absorbers were probably absorbing light 
from background QSOs, making them redder. The effect is possibly 
strong enough to eliminate the QSOs completely from magnitude- 
limited samples, in some cases, since the attenuation in the rest 
frame UV is much stronger than that in the rest frame optical range. 
Akerman et al. (2005) find no evidence for such an effect, but test- 
ing on a larger sample could yield different results. 

The presence of and the nature of dust in QSOALSs is there- 
fore of interest, as it may shed light on the evolution of dust parti- 
cles in cosmic time, may allow the evaluation of selection effects 
and may allow a discrimination between different kinds of dust in 
the intergalactic medium and in the circum-QSO environment. 



1.2 Inference of dust from gas abundances 

Dust has long been suspected to be a constituent of the QSO 
absorbers given the general connections of galaxies with such 
systems, the presence of dust in known galaxies and the evi- 
dence of depletion of gas phase elements onto dust grains in the 
QSOALSs (Meyer, Welty & York 1989; Khare et al. 2004 and 
references therein), similar to that seen in the MW (Hobbs et 
al. 1993; Sembach & Savage 1996; Jenkins 2004), Small Mag- 
ellanic Cloud (SMC, Welty et al. 2001) and Large Magellanic 
Cloud (LMC, Welty et al. 1999a). The signature of dust deple- 
tion is often characterized by the ratio of column densities of 
dominant ionization stages of a refractory (depleted) species to 
that of a volatile (weakly depleted) species. Examples include 
No n/Nzn ii, Nf c n/Nzn ii (Khare et al. 2004 and references 
therein), or N Fc n/Nsi n (Prochaska et al. 2003b), N being the col- 
umn density, in units of number of ions per square centimeter. 
More quantitatively, No n/Nzn n and Nf c n/Nzn n are smaller 
than their cosmic abundance ratios by factors of 2-5 in QSOALSs, 
whereas, in the Galaxy, the values range from 2-100 depending on 
the kind of gas. Akerman et al. (2005) do not find statistically sig- 
nificant difference in Zn metallicity as well as depletions between 
their sample of 15 DLAs in the complete optical and radio absorp- 
tion line systems and the values in optically selected sample pub- 
lished in the literature. 

1.3 Extinction in QSO absorption line systems 

Attempts to quantitatively detect dust by its differential extinction 
of a given object have been difficult. Pei, Fall & Bechtold (1991) 
found higher extinction in a sample of 13 QSOs pre-selected to 
have damped Ly a systems (DLAs), compared to a sample of 15 
QSOs that had no DLAs and concluded that extinction existed in 
the first sample, but not the second: a dust-to-gas ratio of 10% 
of the local interstellar medium was found with the absorber rest 
frame colour excess, hereafter, E(B — V), <0.03. Ellison, Hall & 
Lira (2005) similarly find a slight reddening (E(B - V) < 0.04) 
in 14 out of 42 QSOs with DLAs. The Sloan Digital Sky Survey 
(SDSS) archive has been used for several studies of extinction in 
sightlines to QSOs. Richards et al. (2003) showed that an observer 
frame colour excess, hereafter, A(g — i), which is the difference 
between the actual colours of a QSO and the median colours of 
QSOs at that redshift, could be defined for the SDSS QSO spec- 
tra and could be used to form templates of objects with appar- 
ent degrees of extinction. However, they could not discern if the 
extinction so detected was from the QSO itself or from interven- 
ing systems. Hopkins et al. (2004) used composite QSO spectra to 
show that the extinction towards QSOs is dominated by SMC-like 
extinction, which they argued was predominantly located at QSO 
redshifts. Recently, Wild & Hewett (2005) have found evidence of 
dust in QSOALSs with detected Ca II lines, with E(B - V)=0.06, 
while Wild, Hewitt & Pettini (2005) find strong Ca II absorbers to 
have E(B — V) > 0.1. Murphy & Liske (2004) studied the spec- 
tral indices of QSOs in SDSS Data Release 2 and found no evidence 
for the presence of dust in DLAs at a redshift ~ 3. They derived an 
upper limit of 0.02 on E(B - V). 

1.4 Searches for the 2175 A bump 

The 2175 A bump is a feature of the extinction curve of gas in the 
Local Group. It was discovered in the 1960s using sounding rock- 
ets (Stecher 1965, 1969). The origin of the feature is unknown, but 



Dust and abundances in QSO absorption line systems 3 



it is thought to be made by carbonaceous material in dust or large 
molecules (Draine 2003). It is well correlated with the total density 
of hydrogen (atomic and molecular) and is, by inference, an unsat- 
urated signature of material in dust grains or molecules. Since it is 
known to vary in strength with the steepness of the general extinc- 
tion (Valencic, Clayton & Gordon 2005), knowledge of its behavior 
is an indicator of different sizes and types of grains. Its absence in- 
dicates, in models, a dearth of certain types or sizes of grains in a 
given locale. Understanding its behavior in other galaxies at differ- 
ent times or under different conditions reveals what is generic and 
what is particular about the dust in space. 

Malhotra (1997) reported detection of the 2175 A bump, in the 
extinction curves of QSOALSs, by averaging 92 QSO spectra in the 
absorber rest frames (See section 9.2 of this paper for a discussion 
of this result). Ostriker, Vogeley & York (1990) showed a possible 
detection of the 2175 A bump in a single object. Recently, the pres- 
ence of the bump has been recognized in a few QSOALSs (Motta 
et al. 2002; Wang et al. 2004; Junkkarinen et al. 2004). Yet another 
ISM constituent, the diffuse interstellar bands (Herbig 1995) are 
seen in one system (AO 0235+164, Junkkarinen et al. 2004). 

1.5 Overview of the paper 

In this paper, we proceed to use all three techniques noted above 
to understand dust in QSOALSs: looking for the presence of red- 
dening, searching for the 2175 A bump, and studying the relative 
abundances of elements in a large and homogeneous sample of 
QSOs provided by the SDSS (York et al. 2000). It is of particu- 
lar interest to derive the relationship between these three aspects 
of dust discernment, to test the universality of the characteristics 
of dust throughout the Universe. It is also of interest to distinguish 
dust associated with QSOALSs from the dust near QSOs. Since the 
presence of dust may affect the selection of QSOs for such a study, 
it is important to use samples with well-defined selection criteria 
that can be modeled for the purpose of generalizing conclusions. 
Finally, by surveying well-defined samples of large numbers of ab- 
sorbers, albeit at 170 km s _1 resolution, we can hope to separate 
true trends from selection effects that may arise in small samples 
by accident. 

The contents of the paper, by section, are briefly outlined here. 
Section 2 includes a description of the SDSS survey and spectra, as 
they relate to the present paper. The selection of 809 absorption 
systems for the study of the 2175 A bump is outlined in section 
3. The same sample is used for our derivation of extinction and 
of line strengths. Averaging procedures used to obtain composite 
spectra, along with the line identifications, using the full sample as 
an example, are given in section 4. Two different procedures for 
finding average extinction, one using QSO spectra and the other 
using the SDSS colours are given in section 5. Definition of se- 
lection of sub-samples, based on entities that can be measured in 
individual QSOs, along with the physical motivation of those cri- 
teria appear in section 6. Table 1, found in section 6, summarizes 
the sub-samples and their properties, including the color excesses 
determined as described in section 4. Section 7 contains the cen- 
tral observational results of this study: the effects of extinction are 
correlated with the absorption line strengths and the derived col- 
umn densities; the strength of the Mg II absorption lines is defined 
with regard to extinction and shown not to be a dominant factor; 
and control samples are used to show that there are only small 
effects (on colour excess) in redshift, i-magnitude and beta. The 
extinction is shown to be most like the SMC extinction, with no 
bump and no inflection below 2000 A. The abundances and deple- 



tions for our sub-samples are given in section 8, where we identify 
our highest extinction samples with traditional DLAs (sample of 
~ 100 absorbers) and suggest, assuming an SMC extinction curve, 
that the largest sub-sample (698 absorbers) has, on average, a sub- 
DLA (systems with N H i between 10 19 to 2xl0 20 cm~ 2 ) value 
of Nh i ~ 10 20 cm~ 2 . With these assumptions, the latter sample 
shows solar metallicity in the absorbers. Uncertainties in column 
densities, abundances and extinction values are explored in section 
9. Our conclusions are summarized in section 10. 



2 A LARGE CATALOG OF QSO ABSORBERS 
2.1 The SDSS spectra 

We have undertaken a large survey of QSOALSs, using the QSO 
spectra from the SDSS. The SDSS project uses a 2.5 meter, ded- 
icated telescope at Apache Point Observatory (Gunn et al. 2005). 
The SDSS survey targets QSOs based on colour selection of ob- 
jects to an i magnitude of 19.1 for z em <3.0 and 20.2 for higher 
redshifts, based on 5-filter scans of the northern Galactic cap of the 
MW (Gunn et al. 1998). This choice of sky coverage minimizes the 
effects of Galactic reddening. Based on precise absolute flux cali- 
bration (Hogg et al. 2001; Smith et al. 2002; Ivezic et al. 2004) of 
the object images in the SDSS filter system (Fukugita et al. 1996), 
QSOs can be reliably separated from most types of stars, and rea- 
sonably accurate photometric redshifts can be estimated for a large 
fraction of the quasar candidates, using only colour data (Richards 
et al. 2004; Weinstein et al. 2004). The spectra are obtained using 
fibers plugged into pre-drilled plates, using coordinates known to 
< 0. 1 arcsec (Pier et al 2003). Empirical testing of the technique 
shows that the end-to-end completeness of the QSO spectral sur- 
vey itself is approximately 90%, for reasons discussed by Vanden 
Berk et al. (2005). The completeness is achieved by adaptive tiling 
of selected candidates (Blanton et al. 2003). 

Selection of objects as candidates for taking of SDSS spectra 
occurs in two ways: (a) via a set of rules relating to colour space 
for the purpose of reducing the number of false positives (e.g. white 
dwarfs, A stars) (Richards et al. 2002a) and (b) via the inclusion of 
all known X-ray and radio sources of SDSS objects (Anderson et al. 
2003; Stoughton et al. 2002). Objects selected by the first method 
are brighter than i magnitude 19.1. Objects selected in the second 
way are generally fainter than that magnitude. 

QSOs found using these techniques are entered into QSO cata- 
logs. The catalogs include spectroscopically confirmed QSOs from 
among the candidates described above as well as any other QSOs 
found spectroscopically from among any other targeted objects 
(serendipitous objects, picked as being candidates for some other 
class of point source). The redshifts in the catalogue account for the 
fact that the actual centers of the C IV and Mg II emission lines, for 
example, show systematic shifts with respect to [O III] lines, which 
presumably are at the systemic redshift of the QSO (Gaskell 1982; 
Tytler & Fan, 1992; Richards et al. 2002b). 

The SDSS spectra thus obtained have wavelength coverage 
between 3800 A to 9200 A. The resulting spectra have a resolving 
power of about 1800 (170 km s ) with two-pixel sampling. How- 
ever, the resolution is a complicated function of wavelength, so that 
lines in this paper were observed with a resolution of 80-190 km 
s _1 . The Mg II lines fall in the region of nominal resolution. 

The S/N of the spectra recorded for SDSS targets is set by the 
procedure that 15 minute exposures are recorded until the sum of 
spectra of a single QSO at a g magnitude 20.2 gives (S/N) 2 ~7. 



4 D. G. York, P. Khare et al. 



Spectra admitted to the public database have S/N 4 at a fiber 
magnitude of 19.9 in the red SDSS spectrum and a g magnitude 
20.2 in the blue spectra and a fiber magnitude of 19.9 in the red 
spectrum (Adelman-McCarthy et al. 2005). 

The specific SDSS data products are described by Stoughton 
et al. (2002) and Abazajian et al. (2003, 2004, 2005). The 
data used here are publicly available from the SDSS website, 
http://www.sdss.org. SDSS spectra of QSOs with absorption lines 
have been published in conjunction with broad absorption line 
QSOs (Hall et al. 2002a; Menou et al. 2001; Tolea, Krolik & Tsve- 
tanov 2002; Reichard et al. 2003a) and lensed QSOs (Inada et 
al. 2003; Hall et al. 2002b). All spectra and the photometric data 
used in this paper have been corrected for Galactic reddening, us- 
ing the extinction maps of Schlegel, Finkbeiner, & Davis (1998). 
The QSOs are published in lists, by Data Release number (EDR: 
Schneider et al. 2001; DR1: Schneider et al. 2003; DR3: Schneider 
et al. 2005). Further details of the observing and catalog building 
process may be found in those references. 



2.2 The SDSS absorption line catalog 

Several of us have developed a separate pipeline (Vanden Berk et 
al. 1999; Richards 2001; York et al. 2001; York et al. 2005) to iso- 
late absorption lines and to identify systems of lines. The systems 
are organized into a format modeled on the catalog of York et al. 
(1991). The completeness of the lists of absorption line systems 
is above 95%, as judged by comparing automated lists with those 
made by hand. The incompleteness is associated with the key Mg II 
lines being blended with night sky emission (AA5596, 5890, 5897 
and 6300 and numerous lines above 7000 A), blended with some 
other line of another system or blended together into a single, very 
strong line. Visual inspection was done to verify the present work 
is unaffected by such occurrences. 

The main catalog presents a list of significant absorption 
lines, equivalent widths, various line parameters and possible red- 
shifts. These lines are automatically assigned to discrete systems 
by matching lines with candidate redshifts that agree to within a 
range that depends on the S/N of the spectrum. Weak lines are 
often matched to ±30 km s _1 in the best spectra, whereas strong 
lines in low S/N spectra may be matched to within ±300 km s _1 . 
The systems are graded, depending on the number and quality of 
the individual line detections available. Class A systems are made 
of at least four, unblended lines with detection significance above 
4 a that match in redshift, are free of observational artifacts and 
are not in the Ly a forest. In this paper, we deal only with class A 
systems. Each system has a mean redshift, determined from the av- 
erage of the most significant lines. The redshifts of individual lines 
are retained. Our list of class A systems is complete to the limiting 
equivalent width allowed by each spectrum, for QSOs that meet our 
selection criteria. 

The QSOALS catalog based on DR1 spectra (Abazajian et al. 
2003; Schneider et al. 2003) is used for the current paper, though 
the data are processed with the DR3 software (Abazajian et al. 
2005). 



3 SAMPLE SELECTION 

From a list of ~7,000 grade A absorption line systems extracted 
from 16,713 QSOs in DR1, we used our database to select systems 
with (1) Wm e iia2796, the rest equivalent width of Mg II A2796 > 



0.3 A; (2) 1.0 ^ Zabs ^ 1.86; and (3) an apparent velocity with re- 
spect to the QSO of >3000 km s -1 . The Mg II line strengths were 
used as a first guess to define systems that might have large column 
densities in H I (Rao and Turnshek 2000). The second criterion, the 
redshift range of the absorbers, is to allow only objects for which 
the 2175 A feature lies completely within the SDSS spectrograph 
wavelength range and to ensure that the Mg II lines are at A <8000 
A to avoid the regions of the SDSS spectra that are contaminated 
by strong night sky emission lines. The third criterion rejects asso- 
ciated QSOALSs (Foltz et al. 1986), defined as narrow absorption 
lines within 3000 km s _1 of the QSO redshift. A separate study is 
underway for associated systems. We also imposed an upper limit 
of 1 .95 on the emission redshift of the QSOs to keep the Ly a forest 
from appearing in the spectra. 

These cuts left us with about 1200 systems. For multiple sys- 
tems in the spectra of a single QSO, which are sufficiently close 
in redshift so that their 2175 A features fall within ~ 100 A of one 
another, we only kept the strongest of the systems: this choice max- 
imizes our chances of finding the 2175 A feature in our composite 
spectra without smearing the feature. 

We then examined each of the 1200 spectra by eye. QSOs 
were eliminated if broad absorption line systems (BALs) are 
present; this examination is required to minimize the effects of a 
possible origin of some narrow absorption lines in relativistic out- 
flows from the QSO as well as the possible presence of dust and 
hence reddening in the BALs (Januzzi et al. 1996; Barlow & Sar- 
gent 1997). BALs were identified using the criteria of Weymann et 
al. (1991). In addition, we rejected objects with strong absorption 
systems with relative velocity with respect to the QSO < 3000 km 
s" 1 . 

We were finally left with 813 QSOALSs, out of which we used 
809 systems (for reasons stated in section 5), which we denote as 
sample 1 . A complete list of the 809 systems with relevant data ap- 
pears in the appendix (Tables Ala and Alb). Listed are the QSO 
(plate, fiber and MJD); the coordinates; the QSO redshift; the ab- 
sorption system redshift; i magnitude; j3\ A(g — i) and E(B — V) 
obtained from it (as explained in section 5); and the equivalent 
widths, from the main absorption line catalog, of Mg II AA2796, 
2803, Mg I A2852, Al II A1670, and Fe II A2382. These particu- 
lar lines are among those used in constructing the sub-samples, as 
described later in the paper. 

We assume the spectra are free from contamination of light 
from the host galaxy (Vanden Berk et al. 2005) because we only 
consider QSOs with z em > 1.0. Contamination by the spectrum 
of the host galaxy is well correlated with the presence of an ex- 
tended morphology. We checked that only three of the objects were 
resolved: two appear to be blends with stars and the third, SDSS 
J024634.09-082536.1, is a known, gravitational lens. 

The combination of a nearly fixed signal to noise for each 
magnitude (section 2.1) and the selection of systems with line 
strengths at or above 4 a in significance leaves us with a limiting 
equivalent width that increases with QSO magnitude. We define be- 
low a selection criterion that WM g HA2796 >0.3 A, 4 a. This limit 
is actually reached only in the brightest QSOs. For the faint half 
of our sample, the mean i band magnitude is 19.49; for the bright 
half, it is 18.59. Histograms of the Mg II A2796 line strength show 
that the completeness limit is at an equivalent width of 0.9 A in the 
former case and 0.6 A in the latter case. We will discuss the effect 
of this limit below. 
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Figure 1. The composite spectrum (arithmetic mean) of the systems in the full sample (sample 1) in the absorber rest frame. Individual spectra have been 
normalized using the first 30 eigenspectra of the principle components of the QSO spectra and weighted by inverse variance. All spectra contribute to the 
region between 1900 and 3200 A which, therefore, has high signal-to-noise. Locations of lines of atomic and molecular species are labeled below the spectrum 
while locations of lines of ionic species are labeled above the spectrum. See text for the choice of lines which are labeled. The same figure is shown in the 
Appendix (Figure Al) with a smaller range in flux to highlight the weaker lines. 
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4 GENERATION OF COMPOSITE SPECTRA FOR 
ABSORPTION LINE STUDIES 

Absorption spectra were combined into composites by normaliz- 
ing each spectrum, shifting to the absorber rest frame, masking 
various pixels, weighting, and averaging the spectra. The spectra 
(and associated error arrays) were normalized by reconstructions 
of the QSO continua, using the first 30 QSO eigenspectra derived 
by Yip et al. (2004). The reconstructed QSO continua were formed 
as a linear combination of the eigenspectra in the QSO frame - a 
technique that has the advantage of very closely approximating the 
QSO spectrum while almost entirely ignoring the narrow absorp- 
tion features. After normalizing a spectrum, it was shifted to the 
absorber rest frame, and re-sampled onto a common pixel to wave- 
length scale, which is the same logarithmic scale used by the SDSS, 
but extended to shorter wavelengths. The absorber redshifts used 
were those based on the individual line (e.g. Mg II A2796) that was 
primarily used to select the sub-sample. As noted in section 2.2, 
those may differ from the mean redshift by as much as Az = 0.001 
(300 km s _1 ), but the differences are often as low as 30 kms -1 . 

In several cases, individual pixels were masked before sum- 
ming the spectra. Pixels flagged by the spectroscopic pipeline as 
possibly bad in some way (Stoughton et al. 2002) were masked and 
not used in constructing the composites. Also masked were pix- 
els within 5 A of the expected line positions of detected absorption 
systems unrelated to the target system. The positions of the 18 typi- 
cally strongest transitions were masked, for each system, regardless 
of whether they were actually detected in the spectrum or not. This 
helps prevent strong, unrelated features from "printing through" to 
the final composite, and weaker features from artificially enhanc- 
ing the signal from lines of interest. Such effects would otherwise 
set the noise limits of the composites, especially for the smaller 
sub-samples. 

After masking pixels, the normalized flux density in each re- 
maining pixel was weighted by the inverse of the associated vari- 
ance, and the weighted arithmetic mean of all contributing spectra 
was calculated for each pixel. Error estimates were calculated by 
propagating the weighted errors at each pixel, and by measuring 
the rms scatter about the mean spectrum; both methods resulted in 
similar error estimates. 

Given our choice of redshift range for searching for the 2175 
A feature and including Mg II doublet near 2800 A, the region of 
maximum sensitivity for absorption line detection in the co-added 
spectra is 1900 A to 3200 A in the rest frame. Beyond that range, 
to shorter and longer wavelengths, there are fewer spectra being 
averaged, the noise is larger and the lines (e.g. C IV AA1548, 1550 
and Ca II AA3934, 3969) no longer correspond to the average of 
exactly the same systems as found in the main spectral range. 

Figure 1 shows the region of the composite spectra of 809 
QSOALSs (sample 1) from 1500 A to 4400 A, in the absorber rest 
frame. An expanded version of this figure is shown in Figure Al 
in the appendix A in order to highlight the weak lines. Lines of 
molecules and neutral atoms are labeled below the spectra. Ions are 
labeled above the spectra. We used a complete list of lines in this 
region from Morton (2003) as a reference list. To this we added 
CO lines from Morton and Noreau (1994). The modified list in- 
cludes ample weak lines to confirm detection of each species and 
to find lines of each species that are on the linear portion of the 
curve of growth (lines with strengths oc /A 2 , / being the oscillator 
strength). Only Mg II and Al II, in this spectral region, do not have 
such lines, for our low to medium reddened samples. 

A complete list of lines thus derived is in the appendix (Ta- 



ble A2). As the SDSS absorption line catalog grows, weaker and 
weaker lines will be discernible in the composite spectra. While 
many of the lines are not now detected in QSOALSs, most have 
been detected in the local interstellar medium. The list given here 
provides a minimal set of lines needed to verify identifications in 
these improved spectra. All wavelengths quoted in the paper are 
vacuum wavelengths: when four-digit wavelengths are used in dis- 
cussion, they are truncated values of the wavelengths in Table A2. 

In summary, we are able to co-add many SDSS spectra and 
obtain high signal-to-noise, composite spectra for our selected ob- 
jects between 1900 A and 3200 A in the absorber rest frame. At 
lower S/N, the spectra extend to 1500 A and to 4500 A. The spec- 
tra show interstellar lines from typical, intervening QSOALS. The 
techniques described here are used later to derive spectra of various 
sub-samples of our 809 absorbers. 



5 PROCEDURES FOR DETERMINING EXTINCTION 
5.1 Using composite spectra 

There are numerous ways to define the unreddened continuum for 
derivation of an extinction curve, but division by a QSO spectrum 
with no absorbers is the one we chose, to avoid selection effects due 
to any possible systematic dependence of QSO spectral shape on 
emission redshift or brightness. For the full sample, which we will 
henceforth refer to as the absorber sample, we constructed a match- 
ing sample of QSOs, also referred to as the non-absorber sample, 
having no QSOALSs (of any grade) in their spectra. Thus for every 
QSO in the absorber sample there is a QSO in the non-absorber 
sample which has nearly the same i-band magnitude and nearly 
the same emission redshift (hence, the same absolute magnitude). 
Because the exposure levels are regulated at the telescope, and be- 
cause the exposure levels are checked before a spectrum is admitted 
to the database (see section 2. 1), the matched, non-absorber QSO is 
sensitive to the same level of equivalent widths as the QSO with the 
absorber: systems too weak to be detected in either sample could 
be present in both QSOs of a matched pair, and the effects on ex- 
tinction due to these systems would thus cancel in the derivation of 
the composite extinction curve. 

It is also possible that the QSO spectra are reddened by 
circum-QSO dust. We have not selected the QSOs by any fea- 
ture related to possible extinction other than absorption lines with 
Zabs < z em , so extinction that is circum-QSO should also cancel 
in the composite spectra. 

5.1.1 Selection of the non-absorber sample 

The appropriateness of a QSO with no absorber as a match for a 
QSO with an absorber was determined by minimizing ( ^ Zem ) 2 + 

(%) 2 , R being the rms. Most matches are within Az em < 0.05 
and Ai < 0.1. Gravitational lensing is not expected to significantly 
affect these matchings. Lensing by metal absorbers (Menard 2005) 
and cosmic magnification (Scranton et al 2005) induce, on average, 
magnitude changes that are smaller than our matching range. 

For four of the absorption systems, in three QSOs, it was 
not possible to find a good match (the combined Az em — 
Ai radius is greater than 0.5). The three QSOs are bright 
(i band magnitudes of 17) and are in z em — i regions with 
a low density of QSOs. The objects are SDSS J123602.34- 
033129.9, LBQS0055+0055 (SDSS J005824.75+004113.3) and 
HS 0037+1351 (SDSS J004023.76+140807.3), with emission line 
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redshifts of 1.815, 1.92 and 1.87, respectively, and absorber red- 
shifts of 1.2689, 1.0719 and (in the last object) 1.1254 and 1.6209. 

The sample was thus reduced to 809 systems. We used 174 
non-absorber QSOs twice: none were used three times. The list of 
the matched QSOs along with their i magnitude and z em is given 
in Table Ala (plate, fiber, MJD, emission redshift, i magnitude, 
and match radius) and Table Alb (plate, fiber, MJD, and A(g — 
i)). The quality of the match between the two samples is shown in 
Figure A2, where the difference in i magnitudes for each pair and 
the difference in emission redshifts for each pair are shown. 

Our intent was to form multiple extinction curves using differ- 
ent, non-absorbed QSOs in each case. There are not enough such 
objects in DR1, but it should be possible to use three or four, in- 
dependent, matched pairs with subsequent SDSS data releases for 
absorber samples of the present size. It will also be possible to use 
more sophisticated criteria for matching the QSOs, for instance, 
based on emission line characteristics or the nature of the radio 
flux. 



5.1.2 Derivation of the extinction curves 

Geometric mean spectra were generated from the absorber and 
matching non-absorber spectra. The geometric mean of a set of 
power law spectra preserves the average power law index of the 
spectra, so it is the appropriate statistic to use to determine the 
characteristic extinction law (likely to be approximated by a power 
law) of the absorber sample. The geometric mean composite spec- 
tra were generated in much the same way as the normalized com- 
posites described in section 4, except that the spectra were not nor- 
malized, the pixels were not weighted, and the geometric (rather 
than arithmetic) mean of the spectra was calculated. Weighting the 
pixels can have the effect of changing the final index of the geomet- 
ric mean composite, depending on the way the S/N varies with 
wavelength; we are more interested here in the correct shape of 
the overall continuum than in maximizing the S/N. The geomet- 
ric mean composite spectra were generated for the absorber and 
non-absorber spectra in the same way (both sets of spectra were 
shifted to the rest frames of the detected absorbers). As long as 
the spectroscopic properties of the QSOs in the absorber and non- 
absorber samples are statistically equivalent, except for the effects 
of the absorption systems, dividing the absorber composite by the 
non-absorber composite will yield the relative extinction curve in 
the absorber frame. We use the same procedure to generate extinc- 
tion curves for specific sub-samples described in section 6. 

The geometric mean composite spectra of the 809 absorber 
spectra and the corresponding non-absorber spectra are shown in 
Fig. 2. The composite spectra are smeared out versions of the mean 
spectra as they would appear in the QSO frame. For example, very 
broad features corresponding to the strong Civ, Cm], and Mgll 
emission lines, smoothed and distorted by the distribution of QSO 
vs. absorber redshift differences, are evident in the top panel of 
Fig. 2. Those features are remarkably similar in both the absorber 
and non-absorber spectra, which validates the assumption that the 
two samples are statistically equivalent, at least on scales as large 
as the emission lines. The narrow absorption lines are easily seen 
in the absorber spectrum, and are absent in the non-absorber spec- 
trum. The non-absorber spectrum has been scaled so that the aver- 
age flux density in the reddest 200 A of both spectra are equal. It is 
clear from the figure that the absorber composite spectrum for our 
full sample is relatively fainter at the blue end. 



5.1.3 Nature of the mean extinction curve of the full sample 

The bottom panel of Fig. 2 shows the ratio of the absorber to 
the non-absorber composite spectra. Overlaid are extinction curves 
based on the standard MW curve (Fitzpatrick 1999), and an SMC 
curve (Pei 1992). The SMC curve describes the ratio spectrum quite 
well with a best fit (excluding 5 A around each of the narrow ab- 
sorption features) E(B — l / )=0.013. The statistical error in the 
value of E(B — V) is less than 0.001, and was determined by mea- 
suring the rms dispersion of E(B — V) values of a large sample 
of ratio spectra generated by adding simulated noise to the original 
ratio spectrum. The errors are likely to be larger than this because 
of the fact that the non-absorber sub-samples are not always large 
enough to be completely free of any biases. For example, for the 
sub-samples based on Wing HA2796 values (as described in the next 
section), we expect the A(g — i) values to be the same for all the 
non-absorber samples but they differ by up to 0.03. From the rela- 
tion between < A{g — i) > and E(B — V) obtained in section 9.1, 
we estimate that the errors in the derived E(B — V) values due to 
this in 18 out of 27 samples will be smaller than 0.002. For 7 other 
sub-samples the error will be smaller than 0.004 and only in one 
case (sub-sample 22) it could be as large as 0.009. The MW curve 
also fits the ratio spectrum well using the same value of E(B — V), 
but only at wavelengths longer than about 2500 A. 

There is no evidence of the 2175 A absorption bump in the ra- 
tio spectrum, and the extinction increases at wavelengths shortward 
of 2000 A; both facts are contrary to what is expected for the MW 
extinction curve. We experimented with the generation of compos- 
ite spectra after excising a region up to 8000 km s _1 around the 
major emission lines. We found no effect on the derived E(B — V) 
values. Similarly, using a given matching QSO only once in the 
non-absorber sample by increasing the Az em — At match radius 
also has no effect on the derived E(B — V) values. We note that 
out of all the QSOs in our absorber sample, 65 have two QSOALSs 
while six have three QSOALSs in their spectra. It is possible that 
the reddening is overestimated by the inclusion of these sightlines 
with multiple QSOALSs. To test this we generated composite spec- 
tra excluding these QSOs. The resulting spectrum is almost indis- 
tinguishable (except for somewhat higher noise) from the compos- 
ite of the whole sample. As will be seen below, the systems with 
WMg 11A2796 <L5 A do not, on average produce significant red- 
dening. Out of the 71 sightlines having multiple absorbers, only 13 
have two systems with WMg 11A2796 > 1.5 A. Thus our results are 
not affected by the multiple absorber sightlines. Some of our QSOs, 
picked for having grade A systems, also have grade B (3 detected, 
4 cr lines) and grade C (two lines) or grade D (one line) systems. 
We verified that the grade B systems are subordinate to the grade 
A systems in the strength of Mg II (usually, grade B systems have 
the two Mg II lines and one Fe II line, only) and that the grade C 
systems are mainly C IV-only doublets, with no detectable Mg II. 
The grade D systems are isolated, single lines that could plausibly 
be one member of the Mg II doublet, but the absence of the second 
member makes them insignificant for our study. The effect of the 
presence of systems of lower grades on the results is discussed in 
section 9.1. 



5.2 Using SDSS colours for individual QSOs 

A comparison of the observer frame extinction measure, A(g — i), 
for the absorber and non-absorber samples is shown in Figure 3a 
as histograms. It is clear that the sample with absorption line sys- 
tems as defined herein has more outliers to the right than the non- 




Figure 2. The top panel shows the composite spectrum (geometric mean) 
of the absorbers, in their rest-frames, in red, and that of the corresponding 
non-absorbers, in blue. Note that smeared out emission lines do appear in 
the absorber rest-frame composite spectrum due to the small range in emis- 
sion and absorption redshifts considered here. The absorption lines in the 
red spectra are the obvious, detected absorption lines seen in Figure 1 . The 
bottom panel shows the ratio of the two spectra along with the best fit SMC 
extinction curve in red. The derived absorber rest frame E(B — V) value 
is given. The MW extinction curve for the derived SMC E(B — V) value 
is shown in green. Extinction curves for sub-samples with comparable and 
higher values of E(B — V) are given, in this same format, in Figure 5. 



Figure 3a. Histogram of observer frame colour excess A(g — i) for the full 
absorber sample is shown in red and for the non-absorber sample is shown 
in blue. The fractional number in each bin is plotted. Note the excess of 
QSOs with high values of A(g — i) in the absorber sample. The average 
values of A(g — t) for the absorber and non-absorber samples and full 
DR1 are 0.013, -0.013 and 0.003 respectively. Note that the average for full 
DR1 is close to zero as expected from the definition of A(g — i). The K-S 
probability for the A(g — i) values for our absorber and the non-absorber 
samples to be taken from the same distribution is 1.03xl0 — 5 , while that 
for the absorber sample and the full DR1 is 0.126, and for the non-absorber 
sample and DR1 is 1.15xl0 -6 . The histograms show the clear, red tail of 
our sample of QSO absorbers. 



absorber sample and that the extinction curves we measure mainly 
refer to gas in the intervening absorbers. The nature of the few non- 
absorber QSOs having high values of A(g — i) is the topic of a sep- 
arate paper. These could be due to statistical errors in the colours 
or could involve circum-QSO dust or dust in intervening galaxies 
in an environment with a high dust-to-gas ratio or a high dust-to- 
ri I ratio (i.e., hot gas, with dust, Ferland et al. 2003) or could be 
due to intervening absorbers at z a t s <0.5 (which can not be easily 
detected). They could also involve peculiar energy distributions in- 
trinsic to the QSO itself. Figure 3b shows the same histograms for 
the absorber and non-absorber samples but additionally shows the 
distribution of A(g — i) as a function of z em - The data points are 
colour coded by Mg II line strength. It is clear that the strongest Mg 
II lines can occur in systems of any colour excess, but that the high- 
est relative colour excess objects generally have the highest Mg II 
absorption line strengths as well (the reverse is not true). Out of 1 1 1 
systems in the absorber sample having A(g — i) >0.2, 105 have 

W M g IIA2796 >1.0 A. 

High A(g — i) values appear selectively at high emission red- 
shifts (Figure 3b). This is possible as high z a ts systems are ex- 
pected to be present in the spectra of QSOs with high z em and the 
more reddened (UV) part of the spectrum gets shifted to the SDSS 
filters giving higher values of A(g — i), as A(g — i) is the colour ex- 



cess in the observer frame. We check this as follows. Assuming an 
SMC reddening law, A A = 139\~ 12 E(B - V) (Prevotetal. 1984) 
we determine the rest frame colour excess, E(B— V)( g -i) from the 
observed A(g — i) values (assuming they are due to the intervening 
absorber). Taking X g and Xi to be 4657.98 and 7461.01 A respec- 
tively, this gives E(B - V)( 9 _i) = A(g ~ i)(\+z aha y x/1 1\5§(>. 
In Fig. 3c we have plotted the histograms for E(B — V)( 3 _ 4 ) for 
the absorber and the non-absorber samples. These are similar to the 
histograms for A(g — i) (Fig. 3a) and show that the absorber QSOs 
have higher rest frame E(B — V)( 9 _;) values as compared to the 
non-absorber QSOs. Values of E(B — V)( 9 _;) for individual sys- 
tems, colour coded by the strength of Mg II lines are overplotted as 
a function of emission redshift. As is the case with A(g — i), there 
is still a tendency for higher E(B — with increasing z em . 

This could then be possibly caused by the presence of circum-QSO 
dust in some QSOs as discussed later in the paper. 

Figure 3d shows the same histograms with the absorption line 
redshift plotted instead of the QSO emission redshift. Now, it is 
seen that the distribution of the highest extinction QSOs, in the rest 
frame, is much more uniform. In the last column of Table 1 (below 
in section 6), we give the average E(B — V)r g -i) for the QSOs 
in each of our subsamples. Comparing the color excesses derived 
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Figure 3b. Histogram of observer frame colour excess A(g — i) for the full 
absorber and non-absorber samples are shown in red and blue respectively, 
as in Fig. 3a. Overplotted are the A(cj — i) values as a function of emission 
redshift (labeled to the right) for individual QSOs whose spectra contain 
absorption systems. The points are colour coded to indicate the strength of 
the Mg II A2796 line. Note that the highest values of WM g HA2796 occur 
for all values of A(g — i). 

using our matched pairs technique with the means from the SDSS 
color technique, it is clear that the spectroscopic and photometric 
measures of extinction are very similar. We can use the photometric 
extinction to define sub-samples of objects with similar, high or low 
extinction, which we do in section 6, below. 

In Figure 3e we have plotted the E(B — V)( 9 _;) values as 
a function of Wm s 11x2796- Though the scatter is very large, there 
is a weak correlation between E(B — V)( 9 _ i ) and Wm e iia2796- 
The best fit straight line is E(B - V)=0.016 W Mg HA2796-0.02 
with rms=0.043. The general rise is consistent with the values of 
E(B — V) in Table 1 for the sub-samples defined by ever increasing 
values of Wivig 11A2796, but the scatter of individual points is so 
large that the weak relationship can only be applied to very large 
samples. 

5.3 Summary 

In this section, we have described our matched pair technique for 
the derivation of extinction curves in general and for the full sam- 
ple, in particular; described the selection and uncertainties in the 
matching process; and compared the resulting values with color 
excess derived from SDSS colors for individual objects. The full 
extinction curve for our absorption line selected sample of 809 
QSOALS has no 2175 A bump and no leveling off of extinction 
shortward of the bump. 

We can now proceed to derive spectroscopic extinctions and 
absorption line spectra for sub-samples defined using photometric 



Figure 3c. Histogram of absorber rest frame E(B — V) values for individ- 
ual absorbers as obtained from observed A(g — i) values using the SMC 
extinction curve for the full sample is shown in red. For QSOs in the non- 
absorber sample the absorption redshifts for the corresponding absorber 
QSOs are used to derive the rest-frame E(B — V) values. The histogram for 
the non-absorber QSOs is shown in blue. Note the excess of absorber QSOs 
with high E(B — V)(g_^\ values as compared to the non-absorber QSOs. 
Overplotted are the E(B — V)( g _o as a function of emission redshift (la- 
beled to the right) for individual QSOs whose spectra contain absorption 
systems. The points are colour coded to indicate the strength of the Mg II 
A2796 line. 

measures of extinction in each system, coupled with strengths of 
absorption lines in each system. 



6 SUB-SAMPLE DEFINITIONS 

Using the equivalent widths of the lines discussed above, and other 
parameters, the QSOs were sorted into various subgroups to search 
for any systematic effects in the strengths of lines or extinction de- 
tected in the composite spectra. Table 1 includes the the sample 
number; derived E(B — V) (section 7.1); the selection criteria for 
objects in the sample; the sub-section number describing the sam- 
ple; the number of systems in the sample and the mean values of 
Zabs, /3, i magnitude and E(B — V)( g _j) for the sample. 

For reference, an SMC extinction curve has E(B — V) = 0.01 
for log(N H 1) = 20.6 (Bouchet et al. 1985). 

6.1 Sub-samples in Mg II A2796 equivalent width bins 

The sample of Mg II lines was divided into six ranges of 
WMg 11A2796. ea ch with a similar number of absorbers. Recalling 
that the full sample is called sample 1, we number these six sub- 
samples from 2 to 7: sub-sample 2, 0.3 ^ WMg 11A2796 <0.93; 
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Figure 3d. Same as Figure 3c except that the E(B — V)r g _^ values for 
individual QSOs are plotted as a function of absorption redshift (labeled to 
the right). Note that the distribution of the highest extinction QSOs is much 
more uniform as compared to that in Figure 3d. 

sub-sample 3, 0.93 ^ W Mg 11A2796 <l-2; sub-sample 4, 1.2 < 
W Mg iiA2796 < 1.53; sub-sample 5, 1.53 ^ W Mg 11A2796 <1.91; 
sub-sample 6, 1.91< WMg 11A2796 <2.52; sub-sample 7, 2.52< 
WMg 11A2796 < 5.00 A. Sub-sample 8 has all systems with 
WMg 11A2796 >2.0 A. These groups were defined to explore the 
relation between Mg II strength and the dust content and thereby 
the total gas content if the dust-to-gas ratio is constant. 

6.2 Sub-samples in redshift, magnitude and 8 

The full sample was divided in two equal parts for several param- 
eters: absorbers with z a b s below and above the median value of 
1.3127 (sub-samples 9 and 10, respectively); absorbers in QSOs 
brighter than or fainter than the median value in the i band magni- 
tude 19.12 (sub-samples 11 and 12); and absorbers at /3 below or 
above the median (3, 0.103 (~30,000 km s _1 ) (sub-samples 13 and 
14). These sub-samples were defined, respectively, to look for sys- 
tematic effects in line strength and reddening with redshift, QSO 
brightness, or relative velocity with respect to the QSO. 

6.3 Sub-samples defined to test Routly-Spitzer effect and to 
search for dense clouds 

Similarly, sub-samples were defined by dividing the absorbers into 
two equal sub-sets by line strength: the ratio of equivalent width 
of Fe II A2382 to that of Mg II A2796, with a median value of 
0.577 (sub-samples 15 and 16), and the ratio of equivalent widths 
of Al II A1670 to Mg I A2852, with a median value of 1.538 (sub- 
samples 17 and 18). The first set was designed to see if the Routly- 
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Figure 3e. Plot showing E(B — V)r g _^, the rest frame colour excess as 
obtained from photometric data, for individual absorber QSOs as a function 
of the absorber rest frame equivalent width of the Mg II A2796 line. The 
mean values of E(B — V)r g _^\ with error bars (rms values divided by the 
square roots of the number of points in each of the bins) in Wu g n bins are 
plotted which show the trend of increasing E(B — V)r g —i) with WMg !!• 



Spitzer effect (Routly & Spitzer 1952; Silk & Siluk 1972) is at play, 
wherein the Fe in dust grains might be released in the high velocity 
clouds and strengthen the wings of Fe II relative to Mg II , in a 
manner similar to the strengthening of Ca II with respect to Na I in 
the original effect. The second set was made to look for systems that 
were highly depleted (missing Al II) and of high density (hence, 
Mg II recombined to Mg I). 



6.4 Sub-samples defined by detected, but intrinsically weak, 
lines of heavy elements 

In addition, we defined sub-samples of the systems with detected 
(4 a) lines of the following transitions: Fe II A2260 (sub-sample 
19); Fe II A2374 (sub-sample 20), the blend of Zn II and Mg I 
A2026 (sub-sample 21); and the blend of Zn II and Cr II A2062 
(sub-sample 22). These were defined to test if depletion effects 
could be seen. Fe II might, for instance, become weaker as Zn II 
strengthens, or both might grow, albeit at different rates, with the 
total amount of dust. The Fe II A2374 line is stronger than Fe II 
A2260 line, so yields larger sub-samples to work with. 

6.5 Sub-samples defined to separate Mg II strength from 

A( S - i) 

After examining the results from samples 1 to 22, a set of sub- 
samples was defined based on a combination of Mg II line strengths 
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Table 1. Samples: definitions and properties 



Sample 


E(B - V) 


Selection Criterion" 


Defining 


Number 


< z ahs > 


< P > 


< i > 


< E{B-V) {g _i) > 


number 


(SMC) 




section 


of systems 








1 


0.013 


Full sample 


3 


809 


1.333 


0.120 


19.04 


0.010 


2 


<0.001 


3<Wi, i t <0 93 
* Mg II 


6.1 


129 


1.375 


0.115 


18.66 


0.000 


3 


<0 001 




6.1 


135 


1 318 


1 16 


18 94 


002 


4 


<0 001 


1 ?<W™ TT 53 


6.1 


133 


1 320 


116 


19 06 


-0 004 


5 


007 


1 53<W* T TT ^ 1 Q1 


6.1 


139 


1.332 


0.115 


19.10 


004 


6 


0.018 


1 91<Wu„ tt <2 52 


6.1 


132 


1.316 


0.123 


19.17 


0.021 


7 


0.031 


2.52<W Mg n <5.0 


6.1 


134 


1.327 


0.134 


19.30 


0.035 


8 


0.032 


W U „ TT >2 


6.1,6.5 


251 


1.328 


0.130 


19.26 


0.030 


9 


0.006 


z„ h , <1 3127 


6.2 


404 


1.142 


0.158 


19.02 


0.010 


10 


0.012 


2»k = >1 3127 

^abs ^ * ^ ' 


6.2 


405 


1.523 


0.081 


19.07 


0.010 


1 1 


0.009 


i< 19. 12 


6.2 


398 


1.319 


0.119 


18.59 


0.008 


12 


0.011 


i^l9.12 


6.2 


411 


1.346 


0.121 


19.49 


0.012 


13 


0.010 


/3 <0.103 


6.2 


405 


1.431 


0.052 


19.05 


0.008 


14 


0.011 


/3 >0.103 


6.2 


404 


1.233 


0.187 


19.04 


0.012 


15 


0.009 


W Fe II^Mg II <0.577 


6.3 


368 


1.328 


0.122 


18.96 


0.009 


16 


0.014 


W Fc „AV Mg ii ^0.577 


6.3 


369 


1.299 


0.125 


19.12 


0.013 


17 


0.034 


W A i iiAV Mg i <1.538 


6.3 


85 


1.523 


0.079 


19.04 


0.023 


18 


0.024 


W llll^, IgI ^1.538 


6.3 


85 


1.469 


0.091 


18.89 


0.011 


19 


0.034 


Fe II A2260 present 


6.4 


58 


1.308 


0.126 


19.02 


0.022 


20 


0.019 


Fe II A2374 present 


6.4 


392 


1.300 


0.123 


19.02 


0.014 


21 


0.036 


Zn II-Mg I A2026 present 


6.4 


83 


1.305 


0.132 


18.98 


0.028 


22 


0.058 


Zn H-Cr II A2062 present 


6.4 


31 


1.310 


0.136 


18.76 


0.038 


23 


0.003 


W Mg II <2.0 


6.1,6.5 


558 


1.335 


0.115 


18.95 


0.001 


24 


0.002 


A(g - i) <0.2 


6.5 


698 


1.331 


0.119 


19.02 


-0.002 


25 


0.081 


A(g - i) ^0.2 


6.5 


111 


1.343 


0.125 


19.19 


0.086 


26 


0.012 


W Mg „ ^2.5, A(g -i)< 0.2 


6.5 


97 


1.316 


0.131 


19.34 


0.010 


27 


0.085 


W Mg II ^2.5,A( S - i) ^0.2 


6.5 


48 


1.366 


0.135 


19.27 


0.090 



a For selection from the full sample (sample 1) 
b Equivalent width of Mg II A2796 in A, used throughout in tables 
c Equivalent width of Fe II A2382 in A, used throughout in tables 
d Equivalent width of Al II A 1670 in A, used throughout in tables 
e Equivalent width of Mg I A2852 in A, used throughout in tables 



and A(g — i). We first complemented sub-sample 8 with sub- 
sample 23 (WMg H2796 < 2 A). Then, we divided the full sample 1 
in two parts based on whether the A(g — i) was below or above 0.2 
(sub-samples 24 and 25). Finally, we divided the sub-sample of all 
systems with high Mg II equivalent width (>2.5 A) in two parts, 
according to A(g — i) < 0.2 or A(g — i) > 0.2 (sub-samples 26 
and 27), chosen to sample the extreme reddening values in A(g — i) 
in the histogram of Figure 3a. These sub-samples were defined to 
check if the effects of Mg II equivalent widths and dust could be 
separated, by noticing how the weaker lines in the co-added spec- 
tra depend on the two quantities. 



7 ABSORPTION LINES FOR DIFFERENT SAMPLES 

Average spectra were generated for all 27 samples defined above 
(following the protocols used for Figure 1 (section 4)) and the 
equivalent widths of lines were measured. Then, extinction curves, 
and the corresponding values of E(B — V), were derived for each 
sub-sample, as for Figure 2 (section 5). In section 7.1, we deal with 
the equivalent widths of the lines for a selected set of samples and 
their relation to the derived values of E(B — V). In section 7.2, we 
discuss the relation of the line strengths to other parameters, and in 
section 7.3, we derive column densities and abundances from the 
line strengths. 



7.1 Relation of line strengths to extinction properties 

7.1.1 Separating extinction and line strength effects 

To investigate the systematics that relate extinction to absorption 
line strengths, for our sub-samples, consider the spectra of sub- 
samples 8 and 23-27 plotted in Figure 4. These sub-samples were 
defined to distinguish the effects of large Mg II equivalent width 
from the effects of the SDSS, observer frame reddening, A(g — i). 
The six panels of Figure 4 are arranged in the order of increasing, 
rest frame E(B — V), derived as discussed above in conjunction 
with Figure 2 and as listed in Table 1. 

The absorber, rest frame equivalent widths are listed in Table 
2, for the samples of Figure 4. The header information is the same 
as that on the panels of Figure 4. It is apparent that the Mg II equiv- 
alent widths, while well correlated with the strengths of the strong 
lines of Fe II, do not track well with the spectroscopically derived 
color excesses. The weaker lines, those of Cr II, Mn II and the two 
weakest lines of Fe II, increase monotonically with the color ex- 
cesses. Note the very slow change of the strengths of the corrected 
Zn II equivalent widths (see appendix B) with color excess. 

As is well known, the effects of wide spread, multiple com- 
ponents are strong in Mg II and strong Fe II lines in QSOALS 
(Nestor et al. 2003). However, the effect seems only slightly re- 
lated to extinction. Compare, for instance the two samples with Mg 
II equivalent widths above 2.5 A, one with lower reddening in the 
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Figure 4. Absorber rest frame arithmetic mean composite spectra for sub-samples selected by ranges in A(g — i) and Wmb IIA2796- The panels from top 
to bottom show the spectra for sub-samples 24, 23, 26, 8, 25 and 27, respectively, in order of increasing E(B — V) values. The selection criteria [range of 
Wjv[giiA2796' range of A(g — i)], the number of spectra averaged, the derived value of E(B — V) and the sub-sample number are listed in the lower left of 
each frame. The differing signal-to-noise ratios are due to the number of spectra added together. Only the wavelength of complete overlap in our sample of 
809 systems is plotted; the region below 2300 A is shown on an expanded flux scale to improve the visibility of weak lines in this region. Lines labeled by 
wavelength are as for the same portion of Figure 1. 



absorbers (sample 26, column 6 in Table 2), the other with higher 
levels of extinction (sample 27, column 9). The spectroscopic val- 
ues of E(B — V) show that the extinction differs by a factor 7, 
though both samples have very strong Mg II lines. The extinction 
curves for both samples are SMC extinction curves. If the SMC 
dust to gas ratio applies to both, then the column densities of H I 
are not closely related to the line strengths of Mg II, consistent with 
the findings of Rao, Turnshek & Nestor (2006) (see their Figures 2 
and 7). 

In constructing Table 2, and the later tables, and in anticipa- 
tion of using equivalent widths to derive column densities, we have 
had to face the critical questions of line saturation and component 
blending, especially for the Zn II/Mg I/Cr II blend at A2026 and the 
Zn Il/Cr II blend at A2062. We have used a procedure similar to that 
of Khare et al. (2004) to achieve the separation of these species. The 
detailed process along with the error propagation are given in the 
appendix B, in conjunction with the discussion of column densities 
of our 27 samples. In the process, we have discovered that the Mg 
I A2852 line is unsaturated, for equivalent widths below 600 mA, 
for our samples. 

A second result of the analysis given in the appendix C is that 



the ionization correction is small in all of our samples, as the line 
strengths of Al II and Al III track each other closely, the former line 
implying a larger column density than the latter. 

7.1.2 Extinction and the 2175 A bump 

Following from the discussion in the last section, we can use the 
SDSS photometric colors to isolate the most reddened QSOALS in 
our samples. We now return to our attempt to detect the 2175 A 
bump, which failed for our full sample (Figure 2). Figure 5 shows 
three extinction curves, for sub-samples 26, 8 and 25 of Figure 4, 
for E(B - V) values of 0.012, 0.032 and 0.081. The results are in 
the format of Figure 2. It is remarkable that the typical extinction 
of the MW is so unique compared to our relatively large sample. 
However, we note that virtually all of the known MW extinction 
curves are for relatively large values of E(B — V) (Valencic et al. 
2004) and that for some very low extinction sightlines in the MW, 
the bump is weak or missing (Clayton, Gordon & Wolff 2000). In 
other words, along the very low reddening sightlines between the 
clouds of the MW for which E(B — V) is typically measured, there 
are regions with SMC extinction curves. We have tried to make a 
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Table 2. Equivalent widths for Mg II and relative A(g — i) selected sub-samples 



Sub-sample number 




24 


23 


26 


8 


25 


27 


Number of systems 




698 


558 


97 


251 


111 


48 


Selection Criterion W Ms 


5 II 


all 


<2.0 


>2.5 


>2.0 


all 


>2.5 


Selection Criterion A(g - 


-0 


<0.2 


all 


< 0.2 


all 


>0.2 


> 0.2 




E(B - V) 




0.002 


0.003 


0.012 


0.032 


0.081 


0.085 


A 


Species /A 2 /10 4 




Equivalent width with 1 a errors, in mA 




zojz.yo 


Mgl 


1489 


258±4 


192±5 


596±16 


585±10 


536±15 


795±26 


7H7A AS a 
ZUZO.^+O 


Mgl 


46.4 


<14(>0.5) 


<5 (>8) 


<86 (>24) 


<79 (>20) 


<48 (>21) 


<240(>32) 


1 877 QA 
1 OZ / ,yH 


Mgl 


8.4 


<5 


<5 


<19 


<11 


<24 


<37 


77QA 3^ 
Z / VO. J J 


Mgll 


482 


1362±4 


1076±4 


3035±17 


2674±10 


2198±15 


3210±31 


ZoUj.j J 


Mg II 


241 


1192±4 


927±4 


2712±17 


2400± 1 


2038±15 


2877±30 


10 / u. /y 


A1II 


486 


441±6 


330±6 


1033±24 


1029±14 


943 ±22 


1346±35 


1 77 
1 034. /Z 


Aim 


192 


168±4 


142±4 


332±18 


347±10 


395±18 


622±32 


1 7Q 


Al III 


97 


106±4 


85±4 


198± 16 


221il0 


246± 1 8 


365±33 


loUo.UU 


Sill 6 


0.69 


62±5 


46±5 


187±19 


167±11 


175±24 


212±37 


7fl^£ 7A 


Cr II 


43.6 


35±5 


23±5 


93±16 


101±10 


77±13 


154±27 


2062.24° 


Cr II 


32.4 


23 


17 


65 


76 


75 


131 


2066.16 


Cr II 


21.7 


12±4 


12±4 


28±9 


54±10 


72±18 


109±27 


2576.88 


Mn II 


240 


47±4 


32±4 


143±14 


140±9 


156±17 


237±26 


2594.50 


Mnll 


188 


42±4 


32±4 


108±22 


129±11 


162±19 


277±32 


2606.46 


Mnll 


135 


21±4 


15±3 


60±18 


79±10 


96±18 


147±30 


2382.77 


Fell 


182 


720±4 


540±4 


1786±15 


1600±9 


1359±17 


1975±32 


2600.17 


Fe II 


162 


723±4 


543±4 


1890±17 


1672±10 


1516±16 


2293±30 


2344.21 


Fe II 


62.6 


513±4 


367±4 


1351±14 


1202±9 


1049±16 


1658±29 


2586.65 


Fell 


46 


462±4 


333±4 


1171±15 


1136±9 


1144±17 


1712±29 


2374.46 


Fell 


17.6 


274±4 


195±4 


761±14 


702±9 


624±16 


958±28 


2260.78 


Fell 


1.23 


43±4 


31±3 


137±16 


U3±9 


90±14 


145±18 


2249.88 


Fell 


0.91 


31±4 


19±3 


123±16 


105±10 


113±20 


215±32 


2367.59 


Fell 


0.012 


<9 


<4 


<26 


<16 


<30 


<50 


2012.17 


Coll 


15 


13±4 


9±3 


<26 


56±13 


<30 


<47 


1941.29 


Coll 


12.8 


<4 


<4 


<14 


<8 


<15 


<29 


2026. 14 a 


Znll 


205 


(43) 


(30) 


(52) 


(77) 


(84) 


(57): c 


2062.66" 


Znll 


105 


21 


15 


26 


38 


47 


27: c 


2026.3° 


Zn-Mg 


bl 


35±4 


20±4 


112±14 


Ul±8 


95±15 


266±29 


2062.4 a 


Cr-Zn 


bl 


45±4 


32±5 


91±13 


114±8 


122±14 


158±20 



a Lines so marked are involved in a correction scheme for blends, based on the last two lines of the table, which are the 
measured values. The equivalent widths of the other four lines marked are estimated as described in the text. 



b The Si II A 1808 may be contaminated with S I A 1807 line. However, Jenkins and Tripp (2001) have shown that the ratio of 
column densities of C I and S I is likely to be between 1 6 and 24 for region with Av between zero and 1 .0 in the Galaxy. The 
upper limits and values of equivalent widths of C I lines in our sub-samples indicate that the contamination of Si II A 1808 by 
SI A1807 line is negligible. 

c The Cr II lines are saturated, so the Zn II and Cr II cannot be disentangled due to component saturation. The values marked 
are probably extreme lower limits. See text. 



rough estimate of the fraction of lines of sight that could have MW 
type of extinction curves, by assuming the E(B — V) to be same for 
the SMC and MW lines of sight. We estimate that a maximum of 
30% of the lines of sight could have MW type of extinction curves, 
the bump being masked by the 70% of the SMC type extinction 
curves. 

If we assume that the curves are SMC curves, and that the 
dust-to-gas ratio is the same as in the SMC, then Nh i/E(B — 
V)=(4.4±0.7)xl0 22 cm~ 2 (Bouchet et al. 1985). Then, in order 
of the panels in Figure 4 (top to bottom) and the columns in Ta- 
ble 2 (left to right), the sub-samples have characteristic values of 
Nh i/10 20 (the E(B — V) of each sub-sample is in parentheses) 
of 0.9 (0.002), 1.4 (0.003), 5.4 (0.012), 14 (0.032), 36 (0.081) and 
37 (0.083). According to traditional definition, the last four values 
correspond to DLAs and the first two, to sub-DLAs. 



7.2 Line strength as a function of other parameters 

7.2.7 Line strengths for ranges in Mg II equivalent width 

We have now shown that there is a plausible relationship between 
the E(B — V) of an SMC extinction curve for the absorbers and 
the combination of absorption line strengths and observer frame 
extinction (A(gr — z)). We turn now to samples selected by other 
criteria as out lined in section 6. Table 3 includes the equivalent 
widths for samples with six different ranges of Mg II equivalent 
width, for selected lines of C IV, Mg I, Mg II, Al II, Al III, Si II, 
Mn II and Fe II. For these sub-samples, the strengths of Cr II and 
other, weak lines often give only upper limits for equivalent widths, 
so we have chosen a particular set of lines for which good data 
can be obtained in all six Mg II sub-samples. We include the C IV 
measurements here, as a matter of curiosity, even though we can 
not measure the lines for the lower redshift QSOALSs, meaning 
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Table 3. Equivalent widths for bins in Mg II equivalent widths 



Sub-sample Number 


2 


3 


4 


5 


6 


7 


Number of systems 


129 


135 


133 


139 


132 


134 


Range of Wm 6 


II 


0.3-0.93 


0.93-1.2 


1.2-1.53 


1.53-1.91 


1.91-2.52 


2.52-5.0 




E{B - V) 




<0.001 


<0.001 


<0.001 


0.007 


0.018 


0.031 


A 


Species 


/A 2 /10 4 






Equivalent width in mA 






1548.20 


C IV 


45.5 


928±48 


670±72 


789±64 


1277±48 a 


1 1 08-l-78 a 


1311±90 a 


i W) 78 


C IV 


22.8 


604±43 


543±74 


690±59 


944±51 


SI 8+68 


1 1 S44-1 


2852.96 


Mgl 


1490 


107±10 


171±22 


281±16 


315±26 


491±20 


629±34 


2796.35 


Mg II 


482 


630±19 




1391±20 


1672±25 


99-30-1-99 


1041 +10 


2803.53 


Mg II 


241 


544±18 


ou / in^u 


1167±18 


1455±25 


1979±20 




1670.79 


Al II 


486 


949^40 


98S-M9 


486±45 


664±47 


QflQ+77 


1 09/14-/18 


1854.72 


Al III 


192 


85±16 


1 28± 1 6 


169±15 


209±24 


337±3 1 


407±35 


1862.79 


Al III 


96.5 


43±13 


77±14 


111±15 


136±24 


214±33 


236±27 


1526.71 


Sill 


31 


180±23 


246±35 


319±61 


614±62 


928±73 


1165±126 


1808.00 


Sill 


0.69 


27±13 


63±18 


79±18 


52±16 


147±19 


177±22 


2576.88 


Mn II 


240 


<19 


27±11 


61±14 


58±14 


113±20 


159±19 


2594.50 


Mn II 


188 


<19 


39±14 


43±13 


23±10 


104±20 


130±19 


2606.46 


Mn II 


135 


<19 


<23 


24±11 


<24 


58±16 


93±19 


2382.77 


Fell 


182 


265±10 


5 1 4± 1 9 


714± 14 


910±25 


1330±21 


1728±37 


2600.17 


Fell 


162 


25 1 ± 1 1 


524±22 


775±17 


868±17 


1374±25 


1838±25 


2344.21 


Fell 


62.6 


181±11 


347±18 


521±22 


651±18 




1249±25 


2586.65 


Fell 


46 


138±12 


284±15 


517±21 


564±15 


949±26 


1210±28 


2374.46 


Fell 


17.6 


87±11 


180±18 


281±15 


335±20 


559±17 


679±23 


2260.78 


Fell 


1.23 


25±10 


<22 


58±16 


56±12 


96±16 


176±28 


2249.88 


Fell 


0.91 


<17 


<22 


36±11 


34±12 


55±14 


120±26 



a The strong, blended lines of the C IV doublet had been approximately separated into two lines. The equivalent widths are 
only rough estimates, as the detailed profile is not known. 



that the averaging here is not over exactly the same systems as for 
the Mg II. 

The strong lines of first ions are all saturated (Mg II, Al II, 
one line of Si II, the first three lines of Fe II) and increase in equiv- 
alent width as the mean Mg II equivalent widths increase and as the 
mean extinction increases. Equivalent widths of lines of Al II, Al 
III and Si II A 1526 vary over a wider range (a factor of five as com- 
pared to less than a factor of two for Mg II and C IV), but show a 
monotonic increase with W Mg HA2796- The weak lines vary over 
the largest range (a factor of ten) and show monotonic increase 
with W Mg HA2796 except for a decrease for the sub-sample with 
WMg 11A2796 between 1.53 to 1.91 (sub-sample 5). The weakest 
lines vary by a factor of three while the extinction varies by a factor 
of 30 (compare sub-samples 4 and 7). Evidently, the extinction is 
correlated with the equivalent widths of heavy element lines but in- 
creasing depletion and increasing saturation may preclude finding 
a perfectly linear relationship. 

While we cannot deal with C IV here, as the sub-samples do 
not include the same objects as does the Mg II sub-samples, it is 
interesting to note that C IV behaves differently than the first ions, 
being much less variable across the full range of sub-samples. This 
may be an indication that, in particular, C IV is in a different gas 
than the Al II, and also a different gas than the Al III, since Al II 
and Al III track each other. This result is not inconsistent with the 
literature, but is much more general, because it covers so many ob- 
jects. High resolution data also seems to indicate that Al III follows 
lower ions, not C IV and Si IV. 



7.2.2 Split samples: i magnitude, redshift and (3 

As seen in the last section line strengths appear to increase with 
extinction (Table 2) as well as with Mg II equivalent width (Table 
3). To test whether there is any dependence on other properties of 
QSOALS or QSOs we split the full sample in halves: in absorption 
redshift, at 1.3127; in i-magnitude, at 19.1; and in relative veloc- 
ity with respect to the QSO, at 30,000 km s _1 . The results of the 
line measurements are listed in Table 4, using the same transitions 
that are listed in Table 3. The line strengths are remarkably uni- 
form for samples split in z a bs and /3, implying that they are not 
closely related to the absorption redshift and implying that there is 
no evidence for systems with smaller relative velocity with respect 
to the QSO being different (which could have been the case if they 
were produced in outflows from the QSO) from the systems with 
larger relative velocities (which in all probability are produced by 
intervening absorbers). 

There is a trend for all lines to be stronger in the faint sam- 
ple than in the bright sample. The extinction is the same, however, 
in both samples (11 and 12, see Table 4). Recall that the magni- 
tude that splits the full sample of 809 QSOALS in half is 19.1. As 
mentioned in section 2.1, brighter QSOs are drawn from the SDSS 
color selected sample while fainter QSOs are mainly those chosen 
for spectroscopic follow up because they are X-ray sources, radio 
sources, or serendipitous sources. The bright sample is therefore 
subject to systematic effects of extinction, the faint sample is not. 
There appear, therefore, to be no strong selection effects that hide a 
large number of QSOs due to large extinction, in the SDSS sample. 

The lines are stronger in the faint sample because the weak 
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Table 4. Equivalent widths for split sample: redshift, (3, i magnitude 



Sub-sample Number 


9 


10 a 


U a,b 


U b 


13 a 


14 


Number of systems 


404 


405 


398 


411 


405 


404 


Selection criterion 


2 a u fi <1.31 


2 a i. fi > 1.31 


i<19.1 


i>19.1 


/3 <0.1 


p >0.1 




E{B - V) 




0.006 


0.012 


0.009 


0.011 


0.010 


0.011 


A 


Species 


/A 2 /10 4 






Equivalent width in mA 








C TV 


tJ.J 


NA 


971±23 c 


993±28 c 


948±58 c 


980i27 c 


OJ4-XOO 


1 J JO. / o 


C IV 


22.8 


NA 


727±24 


761±27 


578±39 


732±28 


KJ i J _[_ O J 


2852.96 


Mgl 


1490 


341-j_20 


300±17 


298±13 


379±21 


279± 1 5 


331±12 


2796.35 


Mg II 


482 


1516±14 


1534±14 


1471±13 


1771±17 


1444^14 


1594±13 


2803.53 


Mg II 


241 


1 ^£A4- 1 1. 
1 JDOI LJ 


1 _H)UIL 1 O 


1326±12 


1569±20 


1 252± 14 




1670.79 


Al II 


486 


17 


442±21 


401±26 


522±33 


428±24 




1854.72 


A1III 


192 


190^13 


1 77-f-Q 


164±9 


250±17 


] 79^9 


189il4 


1862.79 


Al III 


96.5 


118±15 


103±10 


95±8 


156±18 


102±9 


137±17 


1526.71 


Si II 


31 


NA 


324±22 


305±25 


397±46 


337±23 


167±55 


1808.00 


Si II 


0.69 


76±14 


58±7 


59±7 


83±24 


59±7 


88±19 


2576.88 


Mn II 


240 


61±8 


53±9 


51±6 


76±17 


43±7 


76±12 


2594.50 


Mn II 


188 


54±7 


42±9 


54±9 


60±12 


47±8 


54±10 


2606.46 


Mn II 


135 


24±7 


<14 


23±8 


<23 


15±7 


29±10 


2382.77 


Fell 


182 


802±10 


710±15 


713±11 


979±20 


695±12 


858±13 


2600.17 


Fe II 


162 


860±11 


742±14 


765±11 


1000±21 


758±13 


884±16 


2344.21 


Fe II 


62.6 


569±10 


488±10 


489±8 


734±17 


486±10 


624±16 


2586.65 


Fe II 


46 


585±11 


477±14 


494±12 


705±23 


494±13 


583±14 


2374.46 


Fe II 


17.6 


307±7 


253±11 


252±8 


416±18 


264±9 


313±9 


2260.78 


Fe II 


1.23 


61±11 


44±10 


48±7 


59±17 


53±10 


42±7 


2249.88 


Fe II 


0.91 


42±8 


19±9 


32±7 


39±12 


30±8 


39±11 



a The numbers of objects included in sub-samples 10, 11 and 13 were larger by 2, 1 and 1, respectively, for the line 
measurements than for the measurements of extinction. The numbers for the number of spectra included and the mean 
extinction are from Table 1 . 

b Note that the lines are stronger in the sample of faint QSOs (sample 12) as the equivalent width limit for detection is higher 
in this sample as compared to the sample of bright QSOs (sample 1 1) 
c Same as footnote a in Table 3. 



systems that are possibly present in that sample are selectively not 
detected at 4 a significance, our cutoff for inclusion. Thus inclusion 
of the faint QSOs may enhance the equivalent widths more than the 
E(B — V) values and therefore may yield somewhat higher col- 
umn densities and abundances. This will affect the column densi- 
ties and abundances measured below, for most other sub-samples 
(which should have an equal distribution of bright and faint QSOs) 
by less than 25% (as seen from the differences in the line strengths 
of samples 1 1 and 12). Also, the trend in the abundance variation as 
a function of E(B — V) that we determine in this paper should be 
free from this effect. Measurements of color excess are evidently 
not affected by the magnitudes of the QSOs. We note, from Table 
3, that samples of equivalent widths of Mg II A2796 < 1.53 A 
(samples 2, 3 and 4) do not have discernible extinction. 

A remaining anomaly from Table 3 is that the derived, spectro- 
scopic extinction for sample 9 (low z abs , E(B — V)=0.006) is one- 
half the extinction for sample 10 (high z abs , E{B - V) =0.012). 
Both extinction values are low. The difference could be significant 
and might be verified by using a larger sample. However, the mean 
extinction is affected by a few high values and many low values. 
Small number statistics in the former is probably responsible for 
the apparent difference. 

In summary, the trends of increasing line strengths with ex- 
tinction (section 7.1.1) and Mg II line strength (section 7.2.1) are 



highly significant compared to only weak trends found by splitting 
the main sample by redshift, (3, or apparent magnitude. 



7.2.3 Other sub-samples 

A final set of sub-samples appears in Table 5, those made to study 
the extremes of the sample from the point of view of the absorption 
lines, as described in sections 6.3 and 6.4. These sub-samples pro- 
vide better limits on the weaker lines and are indicative of a wider 
range of constraints on the QSOALS. 

Table 5 aims to show the weakest lines of each element and 
includes species not represented in Tables 2, 3 and 4: C I, Ca II, 
Ti II, Co II and Ni II. We have also added some strong limits for 
undetected species. 

First, consider Mg I. We have corrected the blends of Mg I, 
Zn II and Cr II as described in appendix B. Here we add the limit 
for Mg I A 1827, which has an / value 20% of that for Mg I A2026. 
The measured limits are consistent with our previous conclusion 
that Mg I is on the linear portion of the curve of growth: for most 
sub-samples we can use A1827 to show that the A2852 line is not 
badly saturated. The case is, however, particularly strengthened for 
sample 1 (the sample of all 809 systems) and for sub-sample 21 
(the Zn II selected sub-sample). The latter has among the strongest 
Mg I A2852 lines of all the sub-samples, 636 mA, a level that for 
any other line in our study would be strongly saturated by a factor 
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Table 5. Equivalent widths for the full sample and for selected sub-samples 





Sample Number 




16 


1 


17 


19 


21 


Number of systems 


369 


809 


85 


58 


83 


Selection criterion 


WthY, TtAV TVTn- TT 
rC Jl "lg 11 


Full 


W A i nAV Mg 1 


Fe II A2260 


Zn II A2026 








>0.58 


sample 


<1.54 


present 


present 




E(B — V) 




0.014 


0.013 


0.034 


0.034 


0.036 


\ 










Equivalent width in mA 




1656.93 


CI 


40.9 


46±13 


19±4 


79±21 


<57 


88±15 


1548.20 


CIV 


45.5 


1 OZXJO 


839±10 


825±44 


970±59 


877-1-47 


1550.78 


CIV 


22.8 


602±37 


580±9 


718±45 


706±62 


551±48 


2852.96 


Mgl 


1490 


7Q74- 1 n 

7 V / XC 1 U 


282±4 


691±35 


689±16 


£7£-L1 ^ 


2026.48" 


Mgl 


46.4 




<32(>11) <37(>26) 


< 74(> 28) 


' i^tu^^z J ) 


1827.94 


Mgl 


8.4 




< 5 


<22 


<15 




2796.35 


Mgll 


482 


1 1 Q_l_ 1 o 
1 O 1 JZE I Z 


1432±4 


1893±28 


2454±16 


7^084- 1 S 

iJUOJI 1 7 


2803.53 


Mgll 


241 


1 UJDT 1 Z 


1261±4 


1751±33 


2255±15 


2312±15 


1670.79 


A1II 


486 




471 ±6 


662±29 


958±26 


936±22 


1854.72 


A1III 


192 


253il3 


183±4 


245±24 


353±16 


345il7 


1862.79 


Aim 


97 


168±14 


112±4 


172±22 


190±15 


202± 17 


1526.71 


Si II 


31 


599±31 


383±10 


510±30 


814±67 


619±55 


1808.00 


Si II 6 


0.69 


121±10 


68±5 


145±22 


240±15 


224±14 


3934.78 


Call 


970 


254±64 


110±16 


NA 


254±47 


J 1 UIJU 


3969.59 


Call 


491 


156±42 


57±13 


NA 


141±36 


154±35 


3384.73 


Till 


410 


67± 16 


31±5 


<130 


160±16 


143±22 


3242.92 


Till 


244 


49Z-IQ 


32±5 


<94 


120±16 


64± 12 


3073.86 


Till 


114 


< 15 


20±6 


<46 


43±12 


<44 


3230.12 


Till 


72 


<? 1 7 


<12 


<61 


<32 




1910.85 


TiII c 


36.5 


^9 


<14 


<19 


32±10 


34±9 


2056.26 


CrII 


43.6 


72± 1 1 


39±4 


89±17 


147±14 


115±12 


2062.24" 


CrII 


32.4 


JO 


40 


70 


1 77 

1Z7 




2066.16 


CrII 


21.7 


39±10 


14±4 


50±16 


100±13 


65±13 


2576.88 


Mnll 


240 


1 14± 1 1 


54±4 


146±21 


23 1 ± 14 


i83±n 


2594.50 


Mnll 


188 


78-1-10 


50±4 


103±16 


167±12 


14.R-M7 


2606.46 


Mnll 


135 


49±9 


26±4 


64±14 


142±13 




2382.77 


Fe II 


182 


1 1 38xx 14 


761±4 


1007±32 


1566±15 


1 S8 1 4- 1 S 


2600.17 


Fell 


162 


1 149± 15 


783±4 


1135±23 


1657±16 


1 VJUV_I_I_ 1 *J 


2344.21 


Fell 


62.6 


M.'JI 1 z 


550±4 


840±22 


1331±14 


1 Z J JZL 1 M- 


2586.65 


Fell 


46 


837±15 


509±4 


845±21 


1447±17 


1 Ton-!- 1 f. 


2374.46 


Fell 


17.6 


505±12 


298±4 


476±31 


894±14 


835±14 


2260.78 


Fell 


1.23 


99±11 


46±4 


125±23 


280±10 


147±12 


2249.88 


Fell 


0.91 


63±10 


35±4 


58±17 


163±15 


110±14 


2367.59 


Fell 


0.012 


<11 


<4 


<17 


<14 


<14 


2012.17 


Coll 


15 


<11 


15±4 


<20 


57±17 


28±8 


1941.29 


Coll 


12.8 


<10 


<8 


< 19 


<27 


< 30 


1741.55 


Nill 


13.0 


84±24 


28±4 


80±20 


124±16 


118±15 


1709.60 


Nill 


9.35 


23±6 


22±5 


<29 


108±22 


65±15 


1751.92 


Nill 


8.58 


<16 


<11 


48±20 


<48 


153±20 


2026.14° 


Zn II 


205 


(42) 


(16) 


(62) 


<148 


<204 


2062.64" 


Zn II 


105 


21 


8 


31 


102 


101 


2026.1° 


Zn II-Mg I 


bl 


55±8 


40±4 


68±17 


176±14 


239±9 


2062.6° 


Cr Il-Zn II 


bl 


77±10 


48±4 


101±15 


225±14 


189±12 



a ' b are same as for Table 2 



c Ti II, A1920.6123, f=0.103 and A1910.9538, f=0.()98, here combined by averaging the wavelengths and adding the f-values. 



of 25. For Mg I, Wivig IA1827 constrains the saturation to be less 
than a factor of 3. 

Next, consider C I. While we have only one line in our spec- 
tral region, it is detected at the 3.5-4 a level in four separate sub- 
samples and is thus real. The abundance of Mg is about 1/10 that 
of C in most known types of interstellar gas (see the later discus- 
sion of abundances). Considering the relative abundances of C and 
Mg and values of /A 2 in column 3 of Table 5, an origin of Mg I 



and C I in the same gas would lead to ratios of the strongest shown 
lines of the two species of three, Mg I A2852 being stronger. The 
observed ratio ranges from 8 to 15 when both lines are detected, 
an amount accounted for by the higher photoionization rate of C 
I compared to Mg I in the interstellar medium. On the other hand 
the recombination coefficient for Mg + compared to C + is much 
higher in warm gas than in gas with T<3000 K (York & Kinahan 
1979). If Mg I and C I are from the same gas, the gas in the four 
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2000 2500 3000 3500 4000 
Absorber Rest Frome Wavelength, A (A) 



Figure 5a. Extinction curve (same format as Fig. 2) for sub-sample 26 (nor- 
malized spectrum which was plotted in panel 3 of Figure 4). The sub-sample 
consists of QSOALS with W Mg IIA2796 > 2.5 A and A(g — i) < 0.2. 
Note that this sub-sample contains only 97 systems as compared to 809 in 
the full sample plotted in Figure 2. Both samples have comparable derived 
E(B - V) values. 

samples that contain C I in Table 5 (16, 1, 17, 21) has a temperature 
lower than 3000 K. Thus, it is consistent that the Mg I and C I arise 
in the same gas. High resolution profiles of each could immedi- 
ately confirm this result, a situation that seems to be true within the 
SMC as well (Welty et al. 2005, in preparation). The importance 
of this conclusion is that the amount of ionized gas in H II regions 
(temperature > 8000 K) is small and the ionization corrections to 
abundances drawn from first ions is small, as is indicated by the 
constant, low ratio of Nai hi to Nai ii- 

To comment briefly on the strong lines, the C IV lines are re- 
markably uniform across these sub-samples, contrary to their varia- 
tion by a factor of two in Table 3. The lines of Al II and Al III show 
a pattern similar to that noticed above; they closely track each other. 
Al II does not track the Fe II A2374 as well as in other sub-samples. 

We turn now to the weak lines and the saturation properties 
of those lines. The last three columns contain the highest average 
values of E(B — V) and there we expect the largest saturation, 
given our analysis of Table 2. However, visual examination of the 
trends and relative strengths in these columns indicates that if we 
take the weakest line in each case, we should get accurate integrated 
column densities. The exception is actually the full sample (sample 
1), the one most dominated by the numerous weak systems, where 
lines of Cr II have equivalent widths ~ 40 mA, instead of the ratios 
expected from the relative /-values of the three lines. However, 
the sizes of the errors for the two, unblended Cr II lines allow the 
possibility that the Cr II lines are unsaturated. Note the apparent 
unsaturated nature of Ca II, in which case the equivalent widths of 
the two lines are similar to the ratio of their /-values, even though 
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2000 3000 4000 

Absorber Rest Frome Wavelength, A (A) 



Figure 5b. Extinction curve (same format as Fig. 2) for sub-sample 8 (nor- 
malized spectrum which was plotted in panel 4 of Figure 4). The sub-sample 
consists of 251 QSOALS with Wj^g UA2796 > 2.5 A (same as for sub- 
sample 26, plotted in Figure 5a) irrespective of their A(g — i) values. The 
inclusion of systems with higher reddening is clearly indicated in the higher 
mean extinction as compared to that for sub-sample 26. 



the total equivalent widths are very large. These lines also have an 
appearance of being broad as do the weakest lines of Mn II and Fe 
II. They may have a wide dispersion of components as in the Mg 
I and Mg II lines discussed in section 9.6, an effect that alleviates 
saturation somewhat. 



7.3 Column densities 

7.3.1 General comments on column density derivations 

We turn now to deriving column densities and abundances from 
the equivalent widths. For the low resolution of the SDSS spec- 
trographs, blends of saturated components are a major concern, 
as discussed earlier and in more detail in the appendix. Jenkins 
(1986) has demonstrated that in large samples, integrated equiva- 
lent widths can be used to derive column densities. Furthermore, if 
multiple lines of a species can be shown to be on the linear portion 
of the curve of growth, the column densities can be directly derived. 
In our best spectra, we can perform this test for lines of Zn II, Cr II, 
Mn II and Fe II. In the extended spectra (1800 A< A and A >3100 
A) lines of Ti II, Co U, Ni II and Ca II are available. Of course, 
when no line is detected, the upper limit on equivalent width can 
reasonably be used to produce upper limits on the column densi- 
ties, provided the data are of high S/N. For the weakest lines, the 
degree of saturation can be assessed if an intrinsically weaker line 
is not detected. 
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Table 6. Column densities 



Absorber somple 
Non-absorber sample 





1.2 

o 

'o 1 

or 



H 1 1 1 1- 



MW Reddening (dashed) 
SMC Reddening (solid) 




2000 3000 4000 

Absorber Rest Frome Wavelength, A (A) 

Figure 5c. Extinction curve (same format as Fig. 2) for sub-sample 25 
(normalized spectrum which was plotted in panel 5 of Figure 4). The sub- 
sample consists of 1 1 1 QSOALS with A(g — i) > 0.2, irrespective of their 
Wjvig HA2796 values. The higher reddening in individual systems leads to 
large derived value of E(B — V). Note the small differences in the strong 
(saturated) absorption lines in panels 3, 4 and 5 of figure 4, but the in- 
creasing strengths of weak lines (towards the left in figure 4) as the mean 
extinction increases. 



7.3.2 Species with weak lines 

We conclude that for our data set, we can derive approximate col- 
umn density averages from Mg I A2852, Cr II A2066, Mn II A2606, 
Fe II A2249 and Zn II A2062 (after correction for Cr II, as in ap- 
pendix B) by assuming these particular lines are on the linear por- 
tion of the curve of growth. Lines of Ti II, Ca II and Ni II can also 
be used in the manner and we assume the same is true for Si II 
A 1808. For cases with evidence of saturation, we treat the derived 
values as lower limits. 

We collect in Table A3, a wide range of sub-samples with 
fairly complete abundance information and with coverage of a full 
range of E(B — V). The table includes equivalent widths for 16 
sub-samples that range in derived E(B — V) from <0.001 to 
0.085. We have purposefully included some sub-samples with sim- 
ilar E(B — V) values, to document the dispersion of equivalent 
widths in differently defined sub-samples. Table A3 includes the 
weakest lines of each detected species to use for abundance esti- 
mates. 





SMC 


Sub-sample 8 


Sub-sample 25 


DLA b 




Skl55 


W Mg n >2.0A 


A(g - i) >0.2 




E{B - V) 




0.032 


0.081 


< 0.04/0 .09 


Species 






Log(N x ) 




HI a 


21.4 


21.2 


21.6 


21.4 


CI 


13.4 


13.1 


13.6: 




Mgn 


16.1 c 


15.0: 


14.9: 




All 




< 12.2 


<12.4 




Al 11 




13.9: 


13.8: 


132 


Aim 


13.0 


13.4 


13.5 


12.8 


Si I 




<12.5 


<12.3 




Sin 


16.2 


15.4 


15.5 


15.6 


CaH 


12.7 


12.4 


12.6 




Ti II 


12.6 


12.4 


12.3 


< 12.9 


CrH 


13.4 


13.3 


13.6 


13.4 


Mn II 


13.0 


12.9 


13.0 




Fel 


<11.6 


<11.8 


<12.0 




Pe 11 


15.2 


15.1 


>14.9 


14.5 


CoH 


< 12.7 


< 13.1 


<13.4 


< 13.0 


NiH 


13.7 


13.9 


>14.1 


14.0 


ZnH 


13.2 


12.6 


>12.7 


12.7 



a For the two sub-samples the H I column densities were estimated from the derived reddening: see section 8 for 
details. For Sk 155 the H I is from a profile lit thai ignores MW gas. from 1 1ST and FUSE data. The column densities of 
heavy elements for Sk 155 refer only to SMC gas. as well. For DLA the H I is from profile fit to Lyman a line from 
HST data. 

b SDSS1727+5302, z abs =1.0311 
c See footnote f to table 7 



Sample W(Mgll 2796) A(g-i) log N 
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▼ Average Upper Limit from Fell 2367 
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Figure 6a. Empirical curve of growth for Fe II for selected sub-samples 
with low (^0.002) E(B — V) values. Column densities are from weak Fe 
II lines (Table A4) and equivalent widths are from Table A3. The curves 
deviate from single component curves of growth because the systems have 
many components of differing column densities. These curves can be used 
to estimate the column densities for systems with only saturated lines avail- 
able (Al II, Mg II, in our case). 



7.3.3 Species with strong lines 

For Al II and Mg II, other techniques must be applied. Figures 6a, 
6b, 6c give curves of growth based on column densities for Fe II for 
the 16 sub-samples of Table A3 with low (Fig. 6a), medium (Fig. 
6b) and high (Fig. 6c) values of E(B — V) respectively. These 
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Figure 6b. Empirical curve of growth for Fe II for selected sub-samples 
with E(B — V) values ~ 0.01. The legend gives the symbols in order of 
decreasing E(B — V), the values of which range from 0.014 to 0.006, 
as listed in Table 1. Column densities and equivalent widths for Fe II are 
derived as for Figure 6a, and the use of the curve of growth is the same. 
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Figure 6c. Empirical curve of growth for Fe II for selected sub-samples 
with high 0.03) E(B — V) values. Column densities and equivalent 
widths for Fe II are derived as for Figure 6a, and the use of the curve of 
growth is the same. 



curves can be used to tie together the weak and strong lines of Fe 
II in an empirical curve of growth. By assuming that the Mg II and 
Al II lines follow the same type of curve (i.e., that the components 
match one to one and have the same relative abundances), the col- 
umn densities of the Al II and Mg II can be estimated. These curves 
of growth have a more linear shape, as saturation is approached 
from the left, than do single component curves of growth. The well 
known effect occurs because lower column density components, 
outside the core absorption, do not affect the integrated strengths 
of the weakest lines, yet add strength to the strong, saturating lines. 
Thus the weak outliers seem to extend the linear portion of the 
curve of growth. We have built such curves of growth for the line of 
sight to SN1987A (Welty et al. 1999a) and to Sk 155 (Welty et al. 
2001), using weak, ultraviolet transitions observed with IUE and 
HST, respectively at high resolution. By adding the column den- 
sity in all components and integrating the entire equivalent widths 
(ignoring the components attributable to the MW), we find that the 
total strengths are similar to those of the DLAs and that the inte- 
grated curve of growth for the sightline is similar to those in Figure 
6. 



7.3.4 An over-view of column densities for several sub-samples 

The column densities along with the estimated errors for the sub- 
samples in Table A3 are given in Table A4. The samples docu- 
mented in this table are chosen to be those that have fairly com- 
plete sets of abundances and to have non-limit values of E(B — V) 
over the widest possible range, 0.002 to 0.085, and to include re- 



dundancy in E(B — V) to measure the scatter from different sub- 
samples. When there is any indication of saturation, we quote the 
formal error based on the assumption that the weakest line of the 
species is on the linear portion of the curve of growth. When sev- 
eral lines of a given species indicate strong saturation, we enter 
values calculated from the weakest lines as lower limits to the col- 
umn density. When there is no detection, we give the upper limit 
to the average column density that could be present. For Al II and 
Mg II the errors are uncertain because of our approximate curve 
of growth corrections. For Si II, we use the values obtained from 
A1808 line only, taken as unsaturated. (The Si IIA1526 line gives 
column densities too low by a factor of 10 on the same assump- 
tion). Errors have been propagated and indicated by a, b, c, d, e,.., 
each representing an increment of 0.05 dex: a implies <0.049 dex, 
b implies <0.099, etc. 

Table 6 includes column densities for two of the sub-samples 
(8 and 25 which are, evidently, typical of the DLA sample in metal- 
licity, depletion and H I column density), taken from Table A4. 
Also given are column densities for Sk 155, the star with the most 
extensive abundance information available in the SMC and for the 
z abs = 1.0311 DLA system toward SDSS 1727+0048 (Khare et al. 
2004). All column densities tabulated are averages over all com- 
ponents, which may differ significantly from one component to 
the other in abundances. For Sk 155, the average is taken after 
the column densities are derived from higher resolution data. For 
the other columns, the data are averages based on lower resolu- 
tion data, as described in this paper. The entries in Table 6 show 
that sub-samples 8 and 25, selected by high Wm e iia2796 and high 
A (g — i) values, respectively, have the same column densities, 
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approximately as the two chosen companion objects, except that 
Nai ii is higher in SDSS 1727+0048 and N Si n is higher in Skl55. 

The column densities for various species from Table A4 are 
plotted as a function of E(B — V) in Figure 7. The general trend of 
increase of column densities with E( B — V) is present but the rel- 
ative increase is small. The increase in column densities is clearly 
not linear in E(B — V), since the single element values range over 
less than a factor of 1 0, while the E ( B — V) values range over a fac- 
tor of 80. It is also clear that the scatter in column density between 
different sub-samples with the same value of E(B — V) is small, 
even though the sub-samples were chosen in very different ways 
and E(B — V) was derived, spectroscopically, after the fact, in all 
cases. We note that an increasing Zn Column density with increas- 
ing extinction has also been noted by Vladilo & Peroux (2005). 

In summary, using equivalent width measurements discussed 
in earlier sections, we elaborate on our technique of finding col- 
umn densities, for elements with multiple weak lines and for Mg II 
and Al II, with strong lines. We find that the column densities fol- 
low a smooth pattern with our E(B — V) values, with little scatter 
between different samples with similar E(B — V) values (Figure 
7). We show that our highest extinction sub-samples have the same 
column density patterns as the SMC and a DLA observed at higher 
resolution. 



8 ABUNDANCES AND DEPLETIONS 

Using the column densities discussed in section 7 (Table A4), and 
the values of H I (also given in Table A4, derived assuming an 
SMC relationship between E(B — V) and H I), we proceed to 
derive abundances, relative to H I, and depletions, relative to zinc. 
The justification for this assumption is that extinction curves for 
all of our samples are of the SMC type. This result may well be a 
consequence of a selection effect associated with the observability 
of QSOALS. There may be other samples of QSOALS for which 
the assumption is not valid. For instance, for the very few systems 
that are known to have the 2175 A bump in their extinction curves, 
the assumption is probably not valid. 

A direct determination of H I for our samples would be prefer- 
able, but as we discuss in section 8.3, this was not possible. We have 
already established on the basis of column densities that our most 
reddened systems can be identified with DLAs at low redshift. Here 
we will show that abundances for sub-samples of smaller extinction 
(sub-DLAs) are close to solar and could only be higher (in case the 
dust-to-gas ratio is higher than the SMC value). 



8.1 Abundances 

The abundances for four sub-samples of Table A4, covering E(B — 
V) from 0.002 to 0.085, derived as explained above, are listed in 
Table 7. The errors quoted are based on propagating the errors on 
the equivalent widths in Table A3. The column densities plotted in 
Figure 7 show that these four systems are a fair sample over this 
range and not special in any particular way. 

In Table 7, it is evident that the abundances, based on Table 
A4 and the SMC dust-to-H I relationship, generally decrease from 
the lowest extinction sub-sample (24) to the highest extinction sub- 
sample (27). The decrease is a factor of 10 to 20. The two elements 
Ca II and C I are trace elements and are not interesting for this par- 
ticular discussion. Comparing the lowest extinction sub-sample 24 
(column 3) with solar system abundances (column 2), we see that 
Zn has cosmic abundances in sub-sample 24. The other elements 
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Figure 7a. Column densities (dex) for Mg II, Mg I, Si II, Ti II, Cr II, Ni II 
and Zn II for all samples from Table A4 as a function of their derived E(B— 
V), assuming the SMC extinction curve applies. Mean values are plotted 
for Mg I. Upper and lower limits are indicated by error bars pointing down 
and up respectively. Upper limits are from non-detections, lower limits are 
from indications of saturation. While the derived reddening differs by a 
factor of 80, the column densities do not increase that much from the left 
to the right. Variations at a given E(B — V) are only a factor of three or 
less. This figure and Figure 7b show that sub-samples defined with different 
selection criterion but having similar values of E(B — V) yield similar 
column densities, a fact we use to justify presenting representative samples 
in Tables 6, 7 and 8. 



are three to ten times lower than the solar system values, Si being 
the least deficient and Mn being the most deficient. On the high 
extinction side, we make an empirical connection to the SMC via 
the star Sk 155 (Welty et al. 2001). The color excesses, H I col- 
umn density and abundances for sample 27 are consistent with the 
values of the Sk 155 and both have an SMC extinction curve. 

Thus, the pattern from our composite spectra is consistent with 
lower total abundances in the more reddened QSOALSs, in refer- 
ence to solar system values, similar to the local, SMC dwarf galaxy 
with sub-solar abundances. Titanium stands out as a possible key 
indicator of different dust depletion patterns that has not yet been 
fully utilized in studies of QSOALSs. It has been detected in a 
few DLAs (Dessauges-Zavadsky, Prochaska & D'Odorico 2002; 
Ledoux, Bergeron & Petitjean 2002). Note that we have not made 
ionization corrections, which are formally required except for Ti II, 
which has the same ionization potential as neutral hydrogen. While 
we argue on the basis of the ratio Nai hi to Nai ii that such cor- 
rections are unimportant, the sense of the corrections would be to 
enhance the relative abundance of Ti noted above, leaving the Ti 
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Figure 7b. Column densities (dex) for C I. Al II, Al III, Ca II, Mn II and Fe 
II for all samples from Table A4 as a function of their derived E(B — V), 
assuming the SMC extinction curve applies. Upper and lower limits are 
indicated by error bars pointing down and up respectively. Upper limits are 
from non-detections, lower limits are from indications of saturation. 



abundances unchanged while lowering the abundances of Zn II, Cr 
II, Ni II, etc. 



8.2 Relative abundances (apparent depletions) 

In Table 8, we compare relative abundances, [X/Zn], a parameter 
most universally available, even when Nh i is not obtainable 1 . The 
solar system values are all zero, in this notation. Also listed are val- 
ues for the "hot", "warm" and "cold" gas in the Milky Way, inferred 
by Jenkins and Savage (1996) and updated by Welty et al. (1999b); 
values for the four samples used in Table 7; values for the SMC 
star Sk 155 (Welty et al. 2001; Welty et al. 2005, in preparation); 
and the average values for five, low redshift DLAs from Khare et 
al. (2004). Nothing in this table depends on the assumption of an 
SMC dust to gas ratio for the sub-samples of this paper. 

The presence of various types of incompleteness: uncertain- 
ties in Mg II and Al II, limits for sample 27 and species that are not 
dominant (C I and Ca II) preclude definite conclusions. However, it 
seems that relative abundances with respect to Zn for sub-samples 
24, 26, and 8 (E(B - V) <~0.032) are most similar to those found 



We note that the comparison should be with Njj=Njj i + Njj 2 + Njj n, 
but the sum is essentially the same as Nh i in all known QSOALSs and 
in the SMC. Ignoring Njj n amounts to our assumption that there are no 
ionization corrections for the dominant first ions. 



in MW sightlines that contain gas with "halo" or "warm" gas de- 
pletion patterns (columns 2 and 3 of Table 8). Only for sub-sample 
27 (E(B - V) ~0.085), Sk 155 and the DLA sample do the rel- 
ative abundances fit into the warm gas category for the MW. In 
none of these sub-samples is there any evidence for "cold" clouds, 
though there well could be "cold components" mixed into the com- 
posite equivalent widths from which we derive the column densi- 
ties. The similarities in these last three columns of Table 8 confirm 
our identification of the most reddened sub-samples with the gen- 
eral class of QSOALSs known as DLAs, as well as our identifi- 
cation of sub-sample 24 with the standard category of sub-DLAs 
(N H i < 2 x 10 20 cm" 2 ). Typical values of H I in the MW halo 
sightlines are N H i ~ 10 19 - 10 20 cm -2 . 

The empirical comparison of relative abundances given in the 
last paragraph does not distinguish the effects of dust depletion 
from those of nucleosynthetic signatures. In regard to the latter, the 
ratio of silicon to iron is the classic pattern found in old Galac- 
tic stars (Lauroesch et al. 1996 and references therein) and the 
value of Mn compared to Fe could be partly due to the odd-even 
effect. However, we do not know that grains are the same every- 
where: Welty et al. (2001) show that the grains in SMC, toward Sk 
155, seem to be mainly Fe and do not contain much Si. The ratio 
[Si/Fe]then mimicks the value found in old halo stars of the Milky 
Way, but two are not related. 

While it is clear from Table 7 that the abundances of first ions 
decrease with increasing E(B — V), the corresponding trend is not 
evident in Table 8, for relative abundances with respect to Zn. In 
particular, the relative abundances of Mg II, Al II, Si II, Cr II, Mn 
II, Fe II and Ni II are lower in the sample with the least extinc- 
tion (sample 24). It will be interesting to look for this trend in high 
resolution studies of large numbers of low z a bs, sub-DLA systems. 

We return now to the issue of our assumption of an SMC ratio 
of E(B — V) to H I column density. If, as might be indicated by the 
relative abundance commented on in the previous paragraph, our 
sub-DLA sample includes some MW extinction curves in the aver- 
aged QSOALS, then there might be less H I per unit of E(B - V). 
In this case, our abundances as given in Table 7 are underestimates. 
The dust to gas ratio can be affected by both grain size and grain 
content. Clayton et al. (2000) show that there are low reddening 
sightlines in the MW that have extinction curves closer to the SMC 
curves as opposed to the high reddening regions of the MW for 
which extinction curves are normally measured. So as long as the 
dust-to-gas ratio in the QSOALS is not lower than the SMC value 
our conclusion about zinc being solar in sub-DLAs appears to be 
robust. 

We conclude that properties that can be known from selected 
QSOALSs from studies of individual QSOs observed at high res- 
olution, are generic to a large sample of QSOALSs with color ex- 
cesses near 0.1. A much larger fraction of QSOALSs than previ- 
ously thought may have solar abundances, judging from the large 
sample size of the lowest extinction sub-sample (number 24) in Ta- 
ble 8. 



8.3 Using known DLA samples for determination of 
dust-to-gas ratio 

We initially considered making use of the H I column densities de- 
rived for a sample of 197 Mg II absorbers by Rao et al. (2005) with 
HST UV spectroscopy, for the analysis of abundances and deple- 
tions. However, we found that their sample did not have significant 
overlap with our high E(B — V) sub-samples. Specifically we di- 
vided their sample into N H i > 10 20 cm -2 and N H i < 10 20 cm" 2 
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Table 8. Comparison depletion in selected sub-samples with other gas 



Gas 






C a 


S24 6 


S26 6 


S8 ft 


S27 b 


SMC C 


DLA d 










-u.z 


-VJ.O 


1 9 


-s. 1 ^ 

S>-l.J 


-V.7 


1 i 


E(B-V) 








0.002 


0.012 


0.032 


0.085 


0.11 




Species 










[X/Zn] 










C 


-0.2 


-0.1 


0.1 


[-0.4] 8 


[-0.6] 


[-0.5] 


[>-0.8] 






Mg 


-0.2 


-0.4 


-0.8 


-0.5: 


-0.2: 


-0.6: 


<-0.5: 


0.0 




Al 




-0.9 


-2.0 


-0.7: 


-0.3 


-0.6: 


<-0.8: 




-0.9 


Si 


-0.2 


-0.3 


-0.9 


-0.3 


+0.1 


-0.1 


<-0.1 


+0.1 


-0.2 


Ca 




-1.8 


-3.2 


K-2.9] 


[>-1.7] 


[-1.9] 


K-1.6] 


-2.2 




Ti 


-0.6 


-1.0 


-2.4 


-0.5 


-0.3 


-0.4 


<-0.7 


-0.9 


-1.0 


Cr 


-0.5 


-1.0 


-1.7 


-0.5 


-0.2 


-0.3 




-0.8 


-0.7 


Mn 


-0.6 


-0.8 


-0.9 


-1.0 


-0.3 


-0.6 


<-0.6 


-1.0 


-0.8 


Fe 


-0.5 


-1.2 


-1.8 


-0.6 


-0.3 


-0.4 




-0.9 


-1.2 


Co 




-0.7 


-1.7 


<0 


<+0.5 


<+0.2 


<+0.4 


<-0.8 




Ni 


-0.5 


-1.2 


-1.8 


-0.6 


-0.2 


<-0.3 




-1.1 


-1.2 


Zn 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



a Representative values for halo, warm and cold Galactic interstellar medium components, from ongoing work noted in the 
text. These values are [X/H], taken here to be [X/Zn] since Zn is solar and undepleted in the local interstellar gas. 
b These columns are the same sub-samples, in the same order, as in Table 7. 
c Values for the star Sk 155, as in Table 7. 

d Depletions for five damped Lyman alpha systems between z abs of 0.8 and 1.5 (Khare et al. 2004). 

e Values in brackets are for species for which there are large ionization corrections, C I and Ca II. Errors for all values are as 
in Table 7. 



Table 7. Abundances™ for selected sub-samples 



Sub-sample Number 


24 


26 


8 


27 




Number of systems 


698 


97 


251 


48 




Criteria 


w Mg II 


All 


>2.5 A 


<2A 


>2.5 A 






A( 9 - <) 


<0.2 


<0.2 


All 


>0.2 




Other 










Sk 155 fa 


E(B 


- V) c 


0.002 


0.012 


0.032 


0.085 


0.11 


Log(N H j) 


20.0 


20.7 


21.2 


21.6 


21.4 


Species d 


Solar e 






Log(N x /N H t ) 






System 












cm 


-3.6 


-7.4 


-7.6 


-8.1 


-8.1 


-8.1 


Mg (II) 


-4.4 


-5.1: 


-5.6: 


-6.2: 


-6.4: 


-5.2^ 


Al (II+III) 


-5.5 


-6.4: 


-6.6: 


-7.3: 


-7.8: 


-7.3: 


Si (II) 


-4.5 


-5.0 


-5.2 


-5.8 


-6.1 


-5.8 


Ca (II) 


-5.7 


-7.8 


>-8.2 


-8.8 


-9.0 


-8.8 


Ti (II) 


-7.1 


-8.1 


-8.2 


-8.7 


-9.3 


-8.7 


Cr (II) 


-6.4 


-7.1 


-7.4 


-7.9 


>-8.0 


-7.9 


Mn (II) 


-6.5 


-7.7 


-7.8 


-8.3 


-8.6 


-8.3 


Fe (II) 


-4.5 


-5.3 


-5.6 


-6.1 


>-6.7 


-6.1 


Co (II) 


-7.1 


<-7.1 


<-7.4 


<-8.1 


<-8.2 


<-8.1 


Ni (II) 


-5.8 


(-6.6) 


(-6.8) 


<-7.3 


>-7.5 


-7.3 


Zn (II) 


-7.4 


-7.6 


-8.2 


-8.6 


>-8.9 


-8.6 



a Abundances arc formally bcllcr limn 0.1-0.2 dcx in precision, assuming ific derived values of E(B — V) can be 
used lo infer Njj j for our samples as described in the text. Values followed by colons are based on saluraled lines 
(Mg 11 and Al If) and Ihe assumption lhal Ihe well determined curve of growth of Fe II can be used to infer N Mg jj 
and jj. Values in parentheses are formally good values but arc uncertain because of inconsistencies with other 

measurements of lines ofNi II. 

k Welty et al. (2001). E(B — V) is the inferred value for the star, corrected for 0.04 mag of extinction due to the 
MW. Nr-r j was derived from a profile fit to the blended Lyman alpha line (from the HST and FUSE data) with SMC 
and MW components. The value given is only for the SMC, as arc the abundances. 
c colour excess for an SMC extinction curve as derived for these sub-samples, as described in the text. 

The ion stage of the listed species from which the abundance is derived is in parentheses. No ionization corrections 
are made, though N^ n jj and N^ n jjj are summed. Ionization corrections are believed not to be important except 
for C I and ca n because N^j jj is greater than N^vj jjj in all of our samples. 
e Solar system abundances from Lodders. K. 2003, Apl, 591, 1220 

f The N Mg jj for Skl55 comes from measurements of the Mg II A A 1238, 1242, with / value 0.0001. It is thus 
much more accurate (formally), but lacking in wings comparable lo the 2800 A doublet. 

sub-samples and made matched non-absorber comparison samples 
using the method described in Section 5. Relative to the comparison 
samples, neither of their sub-samples showed any signature of ex- 
tinction, as both gave E(B — V) slightly less than zero. However, 



this result is consistent with the present work in that it is clearly 
possible to select a sample with a range of Wivig HA2796 that has 
low A(g — 1), as can be seen from Figures 3b and 5a. When Rao 
et al. (2005) selected their sample for HST observation, they chose 
absorbers which had a range of Wm e iia2796 but also required the g 
magnitudes of the QSOs to be brighter than 19.1, thus ensuring that 
there would be sufficient UV flux to obtain the needed signal-to- 
noise ratios to measure Nh i values from the Lya absorption lines. 
By employing this criterion, the most reddened QSOs may have 
been moved out of the their sample. Any small extinction that may 
exist in their sample cannot be easily measured with our technique 
for the small number of objects (15 and 17) in the two sub-samples 
we can check here. 



With regard to the Rao et al. (2005) sample, it should be noted 
that they do find a relationship that can be used to predict the mean 
Nh 1 value as a function of W Mg HA2796- Specifically, they find that 
for W Mg 11A2796 < 0.6A the mean N H 1 ~ 1 x 10 19 cm~ 2 , but for 
Wing HA2796 > 0.6A the mean Nh 1 ~ 4 x 10 20 cm~ 2 . If we were 
to simply use this relationship between Nh i and Wm e 11A2796 we 
would find results similar to those described by Nestor et al. (2003) 
and Turnshek et al. (2005), i.e., a trend of increasing mean metal- 
licity with increasing W Mg 11A2796, where increasing W Mg 11A2796 
is primarily an indication of the increasing velocity spread of the 
absorbing gas. Since neither of these studies separated the effects 
of extinction from those of the kinematics, the result can be thought 
of as a separate effect from depletion, related to our comments in 
section 9.3 on the Routly-Spitzer effect. 
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9 DISCUSSION 

9.1 Extinction and uncertainties 

It is clear that the extinction derived here is due to intervening 
absorbers as it is obtained from the comparison of QSOs with 
absorbers and those without absorbers along their line of sight 
(see the top panels of Figures 2, 5a, 5b and 5c). Additional sup- 
port for this conclusion comes from the correlation of absorber 
properties with the E(B — V) in the following, (i) Extinction 
depends on the strengths of absorption lines e.g. extinction de- 
pends on the mean Wiig HA2796- It is very small for QSOALS 
with W Mg HA2796 < 1 -53 A (sub-samples 2-4) and increases with 
increase in Wm s HA2796 beyond 1.53 A (sub-samples 5, 6, 7 and 
also 8 and 23). (ii) The strengths of the weakest lines, least af- 
fected by saturation, increase directly with the derived values of 
E(B — V) from the extinction curves (see Tables 2, 3, 5 and A3), 
(iii) In sub-samples defined by properties not expected to be related 
to gas absorption lines, the average extinctions and the absorption 
line strengths show little change from sub-sample to sub-sample 
(see Table 4). 

Consistency of the spectroscopic and photometric measure- 
ments of extinction further supports the view that the extinction 
arises in the intervening absorbers in our sample. Both the average 
A(g — i) and average E(B — V")( 9 _j) are correlated tightly to the 
E{B — V) values determined from the composite spectra (the linear 
least squares fits are < E(B - V)( s -o > = 0.98 xE(B - V)- 
0.002 (rms=0.006) and <A( 4 9 '' )> = l.01xE{B - V) - 0.002 
(rms=0.006)). The average value of A(g — i) for the absorber sam- 
ple is positive (0.013) while that for the non-absorber sample is 
negative (-0.013). Also, a higher fraction of QSOs with absorbers 
have large A(g — i) values as compared to the QSOs without ab- 
sorbers. 

To quantify the degree of association of extinction with gas 
absorption, notice that in the absorber sample, roughly 14% of the 
systems have A(g — i) >0.2 as compared to 5% of the systems 
of the non-absorber sample. It thus appears that for 65% of such 
highly reddened QSOs in our sample, intervening absorbers are the 
cause of the reddening. In other words, QSOs with grade A absorp- 
tion systems in their spectra appear to be three times as likely to be 
reddened as those without any absorption system in their spectra. 

As expected the reddening seems to be independent of (3 (sub- 
samples 13 & 14) which is consistent with all the systems with 
(3 > 0.01 being of intervening origin in our BAL-free sample. The 
reddening is also very similar for the sub-samples of bright and 
faint QSOs (sub-samples 11 & 12). Surprisingly, the sub-samples 
selected on the basis of z abs (sub-samples 9 and 10) show very 
different values of E(B — V) \ the sub-sample consisting of lower 
z abs has half the value of E(B — V) as compared to the higher z a b s 
sub-sample. This could in principle, be due to the fact that we have 
only considered systems with grade A in choosing our sample and 
it is possible that there are other systems of lower grade as well as 
systems beyond the range of z abs considered here, along the lines 
of sight. 

However, the only systems not seen by us are those at z abs < 
0.5, which could only be picked up in Ca II. The numbers of these 
systems would be the same in both high z abs and low z abs subsam- 
ples, since all of our QSOs have redshifts above 1. They would also 
be the same, statistically in our matched QSOs, so the effect of such 
systems should not be seen in our extinction curves. As for lower 
grade systems, we have estimated the effect of their presence from 
the redshift dependence of d/V/dz, the number of Mg II systems per 



unit redshift interval, per line of sight, values calculated from the 
SDSS data base (York, Van den Berk et al. 2005, in preparation) 
and find the effect to be < 19%. In reality, the effect is likely to 
be much smaller as our d/V/dz estimates show that the total number 
of systems with grades A, B and C, having W^ig HA2796 > 1.5 A 
(which are likely to cause reddening) is expected to be less than 
< 150 in our sample of 809 QSOs if these QSOs were chosen 
randomly without any knowledge about the presence of absorption 
systems in their spectra. We have, however, chosen these QSOs for 
the presence of absorption systems in their spectra and the probabil- 
ity of their having other absorption systems is thus very small. The 
reason for the large difference between the E(B — V) values for 
the sub-samples with different z a bs ranges is thus not understood. 

9.2 Comparisons with other searches for 2175 A feature and 
dust 

We have established that extinction is from the intervening absorp- 
tion systems and have shown that for none of the sub-samples does 
the extinction curves have a discernible bump at 2175 A, though up 
to 30% of the lines of sight being reddened by MW type of dust can 
not be ruled out. 

There are less than ten known extinction curves for QSOALSs 
that show the 2175 A bump. The best case is AO 0235+164 
(lunkkarinen et al. 2004). Wang et al. (2004) detect the bump in 
three QSOs with z em above 1.6 (out of 6350 in the DR1 sam- 
ple searched). One of these is included in our sample (SDSS 
1144612.98+035154.1). The other two are in QSOs at higher red- 
shift than we include (SDSS 1145907.19+002401.2 and SDSS 
J012147.73+002718.7). The common object has the fourth high- 
est value of A(g — i) in our sample. The three QSOs with the 
highest values of A(g — i) show no sign of the bump, under 
visual inspection (1011117.36+142653.6, 1144008.75+001630.4, 
1085042.24+515911.7). Wang et al. (2004) review previous stud- 
ies, including indications that some lensed QSOs show the bump in 
the spectra of at least one of the multiple images (but, see McGough 
et al. 2005). 

Malhotra (1997) used a subset of spectra from Steidel and Sar- 
gent (1992) to obtain an average extinction curve and to search 
for the 2175 A bump, evidence for which showed up in the com- 
posite spectra. Among the sample are at least two BAL QSOs 
(UM485 (SDSS J121549.8-003422.1) and SBS0932+524 (SDSS 
1093552.92+495314.3)) and several strong, associated absorption 
line systems. Therefore, the sample is not selected in the same way 
as the samples we use here. In addition, some of the spectra have 
strong Fe II emission complexes. 

None of the 19 spectra from Malhotra's sample which have 
been observed by the SDSS, shows the 2175 A absorption feature 
in the SDSS spectra. We did notice that in one of the BAL QSOs 
(UM485), the C IV and Si IV troughs coincide with 2175 A in the 
rest frame of two, narrow Mg II absorption systems. It is possi- 
ble that these BAL troughs contribute to the strength of the feature 
found in the composite spectra. That BAL QSOs may be contribut- 
ing is also indicated by the fact that the composite spectrum ob- 
tained by Malhotra is quite red, a common feature of BAL spectra 
(Reichard et al. 2003b). Thus we believe that the strength of the 
feature found by Malhotra is likely to be caused by artifacts of the 
QSOs chosen, or by a few strong 2175 A features in a few of the 
spectra (ones not yet found in SDSS data). As the SDSS survey is 
completed, a larger number of matches with Malhotra's sample will 
be found and the experiment can be repeated. 

Thus it appears that the carrier of the 2175 A bump is not 
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ubiquitous throughout the Universe, a fact that may go with the ab- 
sence of heavy depletion of heavy elements that we find here in 
our sub-samples. It may be that the bump is not caused by carbona- 
ceous material as is commonly assumed (Draine 2003) but rather 
depends for its existence on heavy depletion of elements such as 
Fe, Ti, Ca and/or other refractory species. 

The ubiquity of the SMC extinction curve also suggests that 
grains in the Galaxy are less completely processed (relatively more 
large grains) than elsewhere in the Universe. Whether the differ- 
ence has to do with selection effects in the QSOALS sample or 
represents global differences between galaxies is not clear. For in- 
stance, dense clouds sampled for MW extinction curves may be 
selectively avoided in QSO spectra and there are low density lines 
of sight in the MW with very weak bumps and very steep UV ex- 
tinction curves (Clayton et al. 2000). 

Ellison et al. (2005) have studied DLAs in a radio selected 
sample of QSOs to overcome any dust bias that may be favoring 
the selection of dust-free absorbers in the magnitude limited QSO 
studies. They obtain a best fit value of E(B - V) = 0.02(0.05) for 
SMC (MW) type dust which is similar, and, in particular, not larger 
than, the values we find here. Their dust-to-gas ratio is however, 
30-50% higher than that for the SMC. 

9.3 The Routly Spitzer effect 

In section 6.3, we defined two sub-samples specifically to test for 
the Routly Spitzer effect (Section 6.3, sub-samples 15 and 16). 
Since the Fe II A2382 line is weaker than the Mg II A2796 line, 
the release of Fe II from grains into gas phase at high velocity 
should yield a high ratio of Fe II to Mg II equivalent widths (us- 
ing the noted lines) if the effect is important. Dividing the sample 
in half (one half with ratios less than 0.58 and the other with ratios 
greater than 0.58) should select for the effect and the equivalent 
widths of other lines may reveal confirming evidence. Comparing 
sub-sample 16, with the higher ratios, with the full sample, there 
is a small effect of enhanced Ca II and Ti II (see Tables A3 and 
A4). However, we have noted elsewhere that several lines seem to 
be partaking in the same wide profiles noted for Mg II. Whether 
the higher velocity components in the broad profiles are selectively 
enhanced cannot be discerned without higher resolution spectra. 

9.4 Interpretation of trends in abundance variations 

Our results show that the gas-phase abundances decrease with in- 
creasing E(B-V) for systems with E(B-V) > 0.002 (the small- 
est measurable value in our sub-samples). Sub-sample 24, having 
698 systems, thus shows the highest abundances. This sub-sample 
has the smallest measurable E(B — V) and hence the smallest 
estimated average H I column density of 10 cm~ and is thus 
composed of sub-DLAs. Peroux et al. (2003) have observed sub- 
DLAs having abundances higher than typical DLAs. Super-solar 
abundance has been observed in two sub-DLAs (Pettini et al. 2000; 
Khare et al. 2004). There are reasons to believe that the DLAs may 
not in fact, be the progenitors of the present day disk galaxies as 
has been suggested earlier (Wolfe et al. 1995) : 

(i) Deep imaging shows that DLAs are often associated with dwarf 
galaxies. 

(ii) The rate of metallicity evolution in DLAs is found to be very 
low and DLAs have been found to have highly sub-solar abun- 
dances even at very low redshifts (Kulkarni et al. 2005a). 

(iii) Star formation is usually associated with molecular clouds, but 



DLAs have very low molecular fraction (Ledoux, Petitjean & Sri- 
anand 2003) 

(iv) The matching of OdlAs at high redshifts with fi star s at present 
may not be correct (Peroux et al. 2005b) or may just be a coin- 
cidence. Also, sub-DLAs may contribute a significant fraction of 
Oh i (Peroux et al. 2005b). 

(v) Zwaan et al. (2005) have shown that the properties and inci- 
dence rates of DLAs are in good agreement with DLAs originating 
in galaxies of various morphological types like the local galaxies. 

On the other hand, Mg II systems, which are likely to be non- 
DLA, Lyman limit systems (LLSs), are very often found to be 
associated with bright spiral galaxies (Bergeron & Boisse 1991; 
Churchill et al. 2005 and references therein). Near-solar abun- 
dances have been observed in one LLS (Jenkins, et al. 2005). Mg 
II systems also show evidence of rotation (Churchill & Vogt 2001). 
Thus it seems likely that non-DLA, Mg II systems, and not the 
DLAs, are the progenitors of present day disk galaxies. The lower 
H I columns in these systems could be a result of higher star forma- 
tion rates as compared to those in DLAs which seem to have very 
low rates of star formation (Kulkarni et al. 2005b). 



9.5 Uncertainties in abundance measurements 

In section 7 and 8 we have discussed several reasons for being cau- 
tious about interpreting the abundances derived from the composite 
spectra. These include the effects of saturation, which in some cases 
cannot be certain and the effects of using the faintest QSOs, which 
have lower equivalent width sensitivity and lead to slight overesti- 
mates of column densities. Other uncertainties are associated with 
the possible presence of circum-QSO extinction without detectable 
absorption lines and with the assumption of an SMC dust to gas 
ratio. Here we discuss these reasons which could affect the abun- 
dance measurements for some sub-samples. 

9.5. 1 The effect of dust around the QSO, with no absorption lines 

First consider the sub-samples based on high A(g — i) values. 
These sub-samples are particularly chosen to contain reddened 
(sub-samples 25 and 27) QSOs. No such criterion is applied for the 
corresponding non-absorber sub-samples. It is of course likely, as 
is clear from the above discussion, that the high values of A(g — i) 
are caused by the presence of absorbers, however, some of the high 
A(g — i) values could have other origins. These sub-samples are, 
therefore, likely to be contaminated (up to 35%) by systems which 
happen to have high A(g — i) values due to reasons unrelated to 
the absorbers. Therefore, even though 65% of the systems in these 
sub-samples may have high dust abundance in the absorbers, the 
remaining 35% may not. Thus, though the E(B — V) for these sub- 
samples will be correct, the equivalent widths will be diluted (as- 
suming for the moment, that high amount of dust in the absorbers 
corresponds to high metallicity and large equivalent widths) due to 
the presence of these stray systems. This will result in getting lower 
calculated abundances for such sub-samples. Since our abundances 
(except for Al II and Mg II) are based on weak lines, this effect is 
of the order of 30% (in the sense of an underestimate) for unsat- 
urated lines, but cannot be evaluated for saturated lines (denoted 
in tables by giving lower limits). For the worst cases, sub-samples 
25 and 27, there are enough reliable determinations that the error 
cannot be large, providing the relative abundances of the heavy el- 
ements do not change drastically at high extinction (see Figure 7). 
The effect is similar in amplitude, but opposite in sign to the effect 
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of including faint QSOs in the composite spectra, as noted at the 
beginning of section 9.5 and in section 7.2.2. 

9.5.2 The assumption of an SMC dust-to-gas ratio 

The largest uncertainties in abundances have to do with the deriva- 
tion of Nh i from E(B — V) using the SMC conversion factor. 
Since the SMC is the extreme curve in all cases known and since 
the less steep extinction curves (such as the standard MW curve) 
have lower values of Nh i for unit E(B — V), the existence of 
some non-SMC curves amidst our sample overestimates the value 
of Nh i. The true, average Nh i would be lower and the resulting 
abundances would be higher than we derive. If anything, our abun- 
dances could be underestimates, though not by too much. 

As argued earlier, there is little evidence for MW extinction 
curves in our sample. Furthermore, Tables 6 and 7 show that sub- 
sample 25 (one of our two most reddened) has column densities and 
abundances (derived by using our technique) that can be related to 
the classic case of damped Lyman alpha absorbers, for which the 
abundances are overwhelmingly sub-solar (see section 1). 

Hence, our conclusion that there is a real difference between 
the low extinction sightlines and our highest extinction sightlines: 
there are solar abundance sightlines (in Zn II) in our lightly red- 
dened, largest sub-samples. The relative depletion of Fe and other 
refractories resembles the halo or warm clouds of the MW, in any 
case, indicating that the Zn effect is due to overall metallicity, not 
to depletion differences on grains in the different sub-samples. 

In the above analysis we have assumed a constant dust-to-gas 
ratio for determining the abundances, and, based on our result of 
finding predominantly SMC extinction curves, used the dust-to-gas 
ratio for the SMC. An examination of all relevant data on the SMC 
reveals little evidence of wide spread variations in the dust-to-gas 
ratio for the SMC (or, for that matter, in the LMC or MW), though 
there are notable exceptions. For 40 known cases in the SMC, for 
instance, the rms scatter in the dust-to-gas ratio is 0.4 dex. If such 
differences existed between a large sample of galaxies with SMC 
curves, they would average out in our composite spectra. It should 
be noted that the stars used in such studies generally have higher 
extinction than our samples with the lowest color excesses and that 
some of the variation for the SMC comes from uncertainties in the 
Galactic foreground extinction. As already noted, low reddening 
sightlines within our Galaxy have curves more similar to the SMC 
curves than to the MW curves derived for higher amounts of ex- 
tinction (Clayton et al. 2000). This variation is consistent with our 
assumptions here. 

In the analysis presented in this paper, we have assumed the 
SMC E(B — V) - Nh i relationship. However, it should be kept in 
mind that if one does not adopt the SMC E(B — V) - Nh i relation, 
but instead takes an approach similar to Nestor et al. (2003) and 
Turnshek et al (2005), i.e., inferring mean Nh i from Win g 11A2796, 
completely different abundance trends are inferred. Of course, a 
variety of trends may be present, representative of different sub- 
populations of absorbers as detailed in section 8 and, below, in sec- 
tion 9.6. Future work with the full SDSS sample may reveal these. 

9.6 Nature of the defining Mg II doublet 

It is believed that large values of WM g 11A2796 are caused by high 
velocity gas. This is also evident from the line profiles in our com- 
posite spectra for various sub-samples. In particular, the strongest 
Mg II lines in sample 7 (2.52 < Wm e 11A2796 < 5.0 A) have average 



widths of 450 km s~ . The correction for the instrumental profile, 
which is typically has FWHM of about 170 km s _1 , is negligible. 
In all probability, some of that width is represented by components 
that are saturated and some by profiles in the wings that are not 
saturated, as seen in the Mg II profile of Peroux et al. (2005a), for 
instance, or as is common in SMC and LMC sightlines. 

The widths of lines of Mg I, Fe II and to a lesser extent, Mn 
II, grow with that of lines of Mg II and it is clear from the line 
profiles that the growth is in added equivalent width over a wider 
velocity range. This fact explains why the absorption lines seldom 
saturate in the sub-samples we measured and why the Mg I A2852 
line is not saturated, despite its large equivalent width, in most of 
the sub-samples measured and presented in Table 2, 3, 4, 5 and A4. 



10 CONCLUSIONS 

We have studied a sample of 809 absorption systems, in the spec- 
tra of SDSS non-BAL QSOs, with Wyg 11A2796 >0.3 A, K 
z abs < 1.86 and [3 >0.01 to determine the average extinction and 
the average strengths of absorption lines in QSOALSs. 

The extinction curves, obtained by comparing the geometric 
mean average spectra of the QSOs in 27 sub-samples of absorbers 
(constructed on the basis of the absorber as well as QSO properties) 
with those of matching samples of QSOs without any absorbers in 
their spectra, show clear evidence of extinction by the intervening 
absorbers. The absorber rest frame colour excess E(B — V) for the 
sub-samples varies from <0.001 to 0.085. The extinction curves 
uniformly share the main properties of the typical SMC extinction 
curve: no 2175 A bump and rising extinction below 2200 A. An 
average MW curve, similar to those of local interstellar clouds, is 
excluded. A mix of up to 30% of MW curves cannot, however, be 
ruled out, but the MW curves are not a dominant contributor. 

Adapting the SMC extinction curve, the observer frame 
colours excesses (A(g — i)) were used, with the absorption red- 
shift, to derive average absorber frame colour excesses, (< E(B — 
V)A(g-z) >), which are consistent with the values of E(B — V) 
derived from the extinction curves. 

While about 14% of the absorber sample has high values of 
A(g — i) (>0.2), 5% of the sample of QSOs with no absorbers also 
shows such high values. These high values, for the non-absorber 
QSOs (and some of the absorber QSOs as well) could be caused by 
strong absorbers with z aba <0.5, where SDSS can not detect Mg 

11 doublets or else, there is a class of QSOs with anomalously red 
continua or there are regions with large amounts of dust associated 
with gas that is too hot to produce the near UV absorption lines. 

The average strengths of absorption lines between 1800 and 
3300 A in the absorber rest frame, are found to increase as the value 
of WMg 11A2796 increases. While Mg II line strengths are generally 
correlated with reddening, there is large scatter: in particular, of 
145 systems with Wjjg 11A2796 above 2.5 A in rest frame equivalent 
width, 2/3 show little reddening (~ 0.01 in E(B — V), on average), 
while the reminder have an average E(B — V) that is eight times 
higher. 

The column densities for Si II, Ca II, Ti II, Cr II, Mn II, Fe 
II, Ni II and Zn II show a strong trend, increasing slowly with our 
derived values of E(B - V). If we adopt the SMC, E(B - V) 
- Nh 1 relation to derive the mean Nh i for our sub-samples, the 
abundances derived from the average column densities and average 
Nh 1 show a consistent trend of lower abundances with increasing 
E(B - V). Our largest sub-sample, with E(B - V) = 0.002 
and Nhi ~ 10 20 cm~ 2 , has near solar zinc and titanium. Sys- 
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terns in this sub-sample (which contains 85% of our full sam- 
ple) are thus sub-DLAs or Lyman-limit systems. For samples with 
E(B — V) ~0.08, we find gas depletions with respect to Zn of up 
to 10 compared to solar abundances. Depletions of Fe, Mn and Cr 
measured with respect to Zn are similar to those for DLAs as pub- 
lished in the literature. The depletion patterns are similar to those 
found in the halo gas in the MW. 

In general, the ionization corrections appear to be small. The 
ratio of the column densities of Al III to Al II is smaller than 0.5. 
The Al III does not seem to be a part of the Al II gas and C IV 
appears to be separated from both Al II and Al III. 

Typical QSOALS have wide blends of velocity components 
in Mg II AA2796, 2802, from below our instrumental width of 170 
km s _1 up to 450 km s _1 or more, independent of extinction. These 
wide systems contribute to the apparent fact that WM g IA2852 is near 
the linear portion of the curve of growth even for equivalent width 
~ 500 mA. Many elements seem to participate in these wide flows 
and the effects of gas depletion seem to be selective, with Ti mildly 
depleted amounts in low reddening samples and highly depleted in 
samples of modest reddening. Since Ti is not yet well observed in 
studies of individual systems, this effect needs to be explored in 
more detail. 

Studies of systems with low but well defined Nh i are needed, 
using high S/N and high resolution, to determine what fraction of 
QSOALS have solar or greater abundances of Zn. 
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Appendix 

A FIGURES AND TABLES 

In this appendix, we collect various figures and tables that support 
the main text, but are illustrative in nature, or not needed in full 
detail in the text but contain important information derived in the 
course of this study. 

Figure Al gives the composite spectrum of the full sample 
of 809 absorbers, but with an expanded relative intensity scale 
compared to Figure 1, to show the weakest absorption lines more 
clearly. 

Figure A2 shows the difference in i magnitude and in QSO 
emission line redshift for the pairs of QSOs matched in the deriva- 
tion of extinction curves. The details are discussed in the text (sec- 
tion 5). 

Table Al contains the list of QSOs used in this study. In both 
Table Ala and Alb, the data on absorber QSOs is on the left and 
that for the matching (unabsorbed) QSO is on the right. Use of 
the pairs is discussed in section 5. Table Ala, includes for QSOs 
with class A absorbers, the plate, fiber and mean Julian date (MJD) 
of the recording of the spectrum for each object; the QSO coordi- 
nates, in decimal degrees; the emission and absorption redshifts; the 
i magnitude; and the relative velocity of the absorber with respect 
to the background QSO. For the matched QSO, the listed items are 
the plate, fiber, MJD, emission redshift, i magnitude and the match 
radius (in units of Figure A2). SDSS data for all objects can be 
found using this information, on the website, http://www.sdss.org 
Any additional information can be obtained from DGY Note that 
the plate and fiber are unique only when the MJD is specified: the 
same QSO can be observed on multiple plates and the same plate 
can be used on different dates. 

Table Alb includes data on the absorbers in the QSOs of Table 
Ala, and on the observer frame reddening of the absorber and non- 
absorber A(g — i), defined in section 1.3. The columns include, for 
the QSO with the absorber; the plate; fiber; MJD; absorption line 
redshift; the value of E(B — V)( 9 _;) inferred from the observer 
frame reddening and redshift, according to the prescription given 
in section 5 (E{B - V) (g _ l) = A(g - i)(l+z a6s )~ 12 /1.506); and 
the rest frame equivalent widths of 5 absorption lines listed in the 
footnotes. For the non-absorber, matched QSO, we list plate, fiber, 
MJD and observer frame reddening. 

Table A2 gives a finding list of absorption lines that may be 
found in spectra such as ours from intervening QSOALSs. Listed 
are the vacuum wavelength, the species, the oscillator strength, and 
the reference from which the data are drawn. The table includes 
only lines between 1370 and 4430 A, and includes all ions of ele- 
ments that can be observed, the main diatomic molecules and the 
main diffuse interstellar band (4430 A) that fall in our spectral re- 
gion. The lines included are chosen to include all lines ever detected 
in QSOALSs as well as additional weak lines that could be used, 
in much better spectra, to confirm single detected lines or to ob- 
tain lines on the linear portion of the curve of growth (e.g., Si II] 
A2335.123). 

The main numerical results of our study are listed in Table 
A3, as the equivalent widths of the weakest absorption lines, useful 
for determining column densities. Sixteen samples of the 27 de- 
scribed in section 6 are included, chosen to cover the full range of 
E (B — V) values we derive and to provide redundancy in sub- 
samples with the same reddening drawn from different selection 
criteria. Four rows across the top include sample number (see Ta- 
ble 1), the number of absorption line systems in each sample, the 
criteria used to select absorbers in the sub-sample and the derived 
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Figure Al. Composite spectra of the full sample; same as figure 1 but the y axis is expanded to show weak lines 
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colour excess for an SMC extinction curve (section 5, Figure 2 and 
Table 1). The vacuum wavelengths, the species and the oscillator 
strengths are given on the left of page 1 of Table A3 and the wave- 
lengths and species are redundantly listed on the left side of page 
two of Table A3, for ease of use. Equivalent widths are given in 
milli Angstroms. Errors are derived from the empirical noise in the 
continuum adjacent to the lines and take into account the width of 
the lines. Since the line widths are comparable to the resolution el- 
ement and the spectra are normalized, using principal components 
for each QSO (section 4), no allowance is made for errors in draw- 
ing the continuum. Errors are 1 a, upper limits are 3 a. Strong lines 
not found in Table A3 may be found, for some systems, in Tables 
2, 3, 4 and 5. 

Table A4 contains the column densities derived using the weak 
lines in Table A3 or, in the case of Al II and Mg II, from the strong 
lines that are available with saturation corrections derived from the 
empirical curves of growth (Figure 6). The column densities tabu- 
lated here are plotted as a function of E(B — V) in Figure 7 of the 
text. 

The banner rows include the same information as listed for Ta- 
ble A3, with the addition of the column density of Nh i, as derived 
from the relation between E(B - V") and H I for the SMC. In the 
second part of the table, the species are redundantly listed down the 
left hand side. 

For each of the sixteen samples, the column densities and up- 
per limits are given in logarithmic form, based in most cases on 
assuming that the weakest line of a species is on the linear portion 
of the curve of growth, for cases where the line ratios justify such 
an assumption. The errors are derived in different ways for dif- 
ferent species. For weak lines and limits, the errors come directly 
from conversion of the equivalent width errors or limits in Table 
A3. The errors are indicated by letters a,b,c etc, each letter indicat- 
ing the logarithmic error, 0-0.049, 0.05-0.099, etc. ForMg II and Al 
II, the additive correction derived from the appropriate Fe II curve 
of growth (Figure 6) is given in brackets. Parentheses, for Ni II, 
indicate that there is some uncertainty beyond that listed because 
different transitions observed do not give consistent values. When 
the assumption of the weakest line being on the linear portion of 
the curve of growth is not justified, the entries are given as lower 
limits (or listed with a colon, if the situation is uncertain). 



B THE SPECIAL CASE OF BLENDED LINES: ZN II, CR 
II AND MG I 

The interpretation of the lines that are blends of lines of separate 
species (mainly Mg I and Zn II near A2026 and Zn II and Cr II 
near A2062) is complicated and we follow procedures similar to 
those outlined by Khare et al. (2004) as follows, (a) Wo 11A2066 
and Wo 11A2056 were averaged, to estimate Wo 11A2062, as the os- 
cillator strength of A2062 is very nearly equal to the average of 
the oscillator strengths of the A2056 and A2066 transitions of Cr 
II. (b) The measured equivalent width for A2062 (last row, Table 
2 as well as Tables 5 and A4) was reduced by this estimate to get 
Wzn HA2062 • This value was multiplied by 2, the ratio of / val- 
ues for the two lines, to get Wz n 11A2026 • The values are in paren- 
theses when the assumption that the Zn II lines are on the linear 
portion of the curve of growth was used, (c) To estimate the con- 
tamination of the Zn II A2026 line by Mg I, we calculated the value 
for W Mg 1A2026 from the linear curve of growth assumption, using 
Wm s ia2026=Wm s IA2S52/42 1 . We compared this with the residual 
left for Mg I after subtracting the Zn II A2026 estimate. We found 
that the calculated contribution of Mg I is remarkably close to the 
lower limit assuming the linear part of the curve of growth applies 
for Zn II A2026. This interesting result, discussed later, is consistent 
with the failure to detect Mg I A 1827 which would be detectable if 
Mg I A2026 were highly saturated. To be conservative (see dis- 
cussion below about blends of saturated lines), we actually enter 
the upper limit to Mg I A2026 by showing, on the left for Mg I 
A2026 in Table 2, the maximum value of Wm s IA2026 assuming 
that Wzn iiA2026=Wzn 11A2062, that is, that the zinc is fully satu- 
rated. The entry to the right (in parenthesis) for Mg I is the mini- 
mum value given the A2852 equivalent width. Only for sub-sample 
27 is Mg I the dominant contribution to the strength of the A2026 
line; this sub-sample is the only one with the equivalent width of 
A2852 > 600 mA, in the table. 

Any conclusions on abundances are strongly affected by our 
technique to unscramble the Zn Il/Cr II/Mg I equivalent widths. If 
the Cr II A2062 line is saturated on a component-by component ba- 
sis, so probably is the Zn II A2062 line. At any rate, component-by- 
component saturation of one species in a complicated blend means 
there is no way to separately extract information on the second 
species. Thus, in the last two sub-samples, the equivalent widths of 
Zn II lines do not grow with E(B — V), possibly because we have 
assigned most of the equivalent width to Cr II using our procedure 
and Zn II is thereby, underestimated. This problem is accentuated 
by the large probable contribution of the Mg I A2026 line, where 
neither Wz„ 11A2026, as we predicted by using our simple rule, nor 
the linear curve of growth prediction for W Mg 1A2026 comes close 
to explaining the 266 mA measured equivalent width of the blend: 
Mg I must be saturated, consistent with the much larger value of 
WMg IA2852 compared to the rest of our sub-samples: in this case 
the Zn II/Mg I blend can not be separated either. 



The errors propagated in these steps are as follows: (a) The error in 
WoiiA2062 = 1A/2 X the error in W CrII A2056; ( b ) The error in 
W Zn llA2062 = X the error in W011A2062; (c) The error in 

W Zn IIA2026 = 2 X the error in W Cr HA2056- The error in W Mg IA2026 
should be very small, if the Mg I A2852 line is on the linear portion of the 
curve growth 
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Figure A2. Plot between difference in z e m and i magnitude of the absorber QSOs and corresponding non-absorber QSOs. Points marked in red are excluded 
from the sample as the differences are large. 



C USE OF BOTH WEAK AND SATURATED LINES IN 
DERIVING COLUMN DENSITIES 

It follows from Appendix B that both Zn II lines are near the linear 
portion of the curve of growth in a number of cases and that we 
can derive Nz n n from the corrected value for Wz„ ha2062- This 
assumption should be verified by getting high resolution observa- 
tions of the four lines involved here for a number of sources in our 
sub-samples. 

Turning to other weak lines, Cr II is consistent with no sat- 
uration in the first four sub-samples, through sub-sample 8, with 
E(B — l/)=0.032, but strong saturation is evident in the sub- 
samples in the last two columns (E(B — V) ~0.08). Mn II is 
consistent with little saturation for the first three sub-samples, but 
shows saturation in the sub-samples in the last two columns. The 
weak lines of Fe II show similar effects. Note that when the weakest 
lines of Fe II, the lines of Cr II, or the Mn II triplet become satu- 
rated (sub-samples 25 and 27), the large errors in equivalent width 
for the sub-samples allow values that seem to be inverted. 

Finally, a pattern related to saturation is discernible assuming 
the first ions lie roughly on the same curve of growth in a given 
average spectrum (Spitzer & Jenkins 1976). Inspection of Table 
2 shows that lines can be assumed to be unsaturated when their 



equivalent width is <50 mA. Lines of higher equivalent widths may 
or may not be saturated. 

We note the close tracking of Al II A 1670 and Al III A 1854. 
This implies that the average ionization correction does not vary 
much between the sub-samples as defined, though the two ions may 
well be in different physical regions, proximate to each other but 
not mixed. The fact that Al II A1670 tracks Fe II A2586 shows that 
it may be possible to estimate Al II column densities by using the 
Fe II curve of growth, which is well defined because of the large 
number of Fe II lines. Interestingly, Al III shows modest saturation 
at all apparent values of E(B — V), consistent with its being in 
separate regions that do not contribute much to reddening but do 
have appreciable column densities. Thus, ratios of Nai iii/Nai ii 
are upper limits to the ratios in the H I regions. Our samples indi- 
cate Nai iii/Nai ii < 0.5 (see below), consistent with the findings 
of Howk & Sembach (1999) and Vladilo et al. (2001). Our large 
sample makes this result very general. 
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Table Ala: Full sample of the absorber and non-absorber QSOs 
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Table Alb: Full sample of the absorber and non-absorber QSOs 
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1 0177 
1 .Uj / / 


O OSSA 
U.UJ J^f 


O 01 S7 
U.U 1 J / 


79 
U. /Z 


1 SI 


1 70 

3. I\) 


u 


1 17 
1.3/ 


OAOA 


SSA 


S 1 8 1 9 
J lo 1Z 


1 191 
-U. 1 JZJ 


097ft 

UZ 1 U 


91 1 
Z 1 J 


S 1 QOQ 


1 1970 


9988 
U.ZZoO 


OSS1 


98 
U.Zo 


1 QQ 

1.77 


1 7S 


QS 


1 18 
1 - JO 


09 Q6 
uzyu 


A17 
H-l / 


S1Q8A 
J Lyo'-t 


1 1A6 
-U. 1 JH-U 


09 7ft 

UZ 1 U 


907 

Z7 / 


S 1 QOQ 
J 1 y\)y 


1 S601 


9869 
u.zooz 


06 1 S 


2 21 


6 60 


6 SS 


2 18 


9 SQ 

Z.J7 


OSAQ 


918 

ZJo 


S1Q81 
j I701 


1790 
-u. 1 / zu 


0970 
UZ 1 U 


462 


S 1 QOQ 
J 1 y\)y 


1 96AQ 


0819 


0907 


40 


1 20 


1 08 


n 
u 


78 
u. / 


OSS9 

UJ JZ 


1 SI 
1 j j 


S 1 QQ9 

j 1 yyz. 


9171 
-u.z 1 / J 


0271 


1 SA 


S1 881 

J 1 OOJ 


1 2240 


0701 

U.U / Uj 


017Q 
u.ui ly 


SO 


1 61 
1.01 


1 11 


n 
u 


8A 
U.ot 


OA 10 


072 


S 1 877 

j 1 / / 


0808 

-U.UoUO 


0273 


1 Q7 

Ly 1 


5 1957 


1 0760 


2709 


0749 


91 


2 46 


9 1 1 

Z. 1 3 


n 
u 


1 42 


0283 


560 


S 1 QSQ 

j 1 yjy 


OOOQ 


0971 
UZ / J 


AQA 


S 1 QS7 

j 1 vj / 


1 0170 


1 01 ft 

-U. IUjO 


09 Q 6 

-u.uzyo 


U 


1 Q9 

1 .yz 


1 S6 

1 .JO 


u 


77 

u. / / 


01S1 

U J J 1 


ASO 
toU 


S 1 780 
J 1 / ou 


O 0616 
-U.UOJO 


0273 


496 


51957 


1.1936 


0.0883 


0.0228 





0.70 


0.97 





0.62 


0650 


185 


52143 


-0.0339 


0274 


357 


51913 


1.4243 


0.0822 


0.0189 


0.38 


2.52 


2.10 


0.63 


1.41 


0533 


029 


51994 


-0.2390 


0275 


357 


51910 


1.0112 


0.0209 


0.0060 





0.65 


0.72 





0.44 


0552 


232 


51992 


-0.0348 


0276 


068 


51909 


1.2846 


-0.0697 


-0.0172 


0.19 


0.84 


0.83 


0.43 


0.46 


0496 


338 


51988 


-0.1248 


0276 


579 


51909 


1.0785 


0.0784 


0.0216 





1.14 


0.97 





0.39 


0315 


598 


51663 


-0.1361 


0276 


579 


51909 


1.7858 


0.0784 


0.0152 


0.27 


1.87 


1.70 


0.73 


0.99 


0315 


598 


51663 


-0.1361 


0277 


059 


51908 


1.6405 


-0.0901 


-0.0187 





0.58 


0.80 





0.40 


0484 


119 


51907 


0.0623 


0277 


246 


51908 


1.1623 


0.0283 


0.0074 


0.75 


1.31 


1.18 





0.74 


0463 


076 


51908 


0.0301 


0277 


417 


51908 


1.6811 


0.1806 


0.0367 





1.02 


0.91 





0.47 


0657 


161 


52177 


0.1214 


0277 


628 


51908 


1.3352 


-0.0000 


0.0000 





1.49 


0.74 





0.58 


0382 


024 


51816 


0.0858 


0277 


629 


51908 


1.0198 


-0.1291 


-0.0369 





1.08 


0.77 





0.70 


0282 


488 


51658 


-0.0653 


0278 


009 


51900 


1.0134 


-0.0521 


-0.0149 





0.88 


0.78 





0.56 


0539 


484 


52017 


-0.0836 


0278 


446 


51900 


1.3274 


-0.1470 


-0.0354 


0.22 


1.34 


1.01 


0.50 


0.52 


0513 


125 


51989 


0.0555 



a Equivalent width of Mg I A2583 A 
6 Equivalent width of Mg II A2796 A 
c Equivalent width of Mg II A2803 A 
d Equivalent width of Al II A 1670 A 
e Equivalent width of Fe II A2382 A 
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Table A2: Finding list for interstellar lines in spectra of QSOs: 1370-4500 A 



Vacuum Species Oscillator Reference Vacuum Species Oscillator Reference 

Wavelength strength Wavelength strength 



1370.132 


Nill 


7.690E-2 


M03 


2325.4029 


CII] 


4.780e-8 


M03 


1381.476 


PI 


3.160E-1 


M03 


2326.1126 


CII]* 


5.520e-8 


M03 


1393.324 


Nill 


1.010E-2 


M03 


2328.8374 


CII]* 


2.720e-8 


M03 


1393.7602 


Si IV 


5.130e-l 


M03 


2335.123 


Sill] 


4.250e-6 


M03 


1402.7729 


Si IV 


2.540e-l 


M03 


2344.2139 


Fell 


1.140e-l 


M03 


1477.222 


Nill 


9.720e-4 


M03 


2346.2616 


Nil 


1.430e-l 


M03 


1477.57 


CO APXS2 


4.013e-2 


MN 


2367.5905 


Fell 


0.216e-4 


M03 


1480.9546 


Co II 


1.190e-2 


M03 


2367.7750 


All 


1.060e-l 


M03 


1509.75 


COAPXS1 


2.91 le-2 


MN 


2374.4612 


Fell 


3.130e-2 


M03 


1526.7070 


Si II 


1.330e-l 


M03 


2382.7652 


Fell 


3.200e-l 


M03 


1532.5330 


PII 


3.030e-3 


M03 


2407.9875 


Col 


3.800e-l 


M03 


1533.4316 


Si II* 


1.320e-l 


M03 


2425.6708 


Col 


2.820e-l 


M03 


1544.45 


CO APXS0 


1.578e-2 


MN 


2448.4509 


Fel 


2.310e-3 


M03 


1547.9451 


Co II 


1.110e-2 


M03 


2463.3922 


Fel 


5.320e-2 


M03 


1548.204 


CIV 


1.899e-l 


M03 


2484.0209 


Fel 


5.440e-l 


M03 


1550.781 


CIV 


9.475e-2 


M03 


2501.8858 


Fel 


4.930e-2 


M03 


1552.7624 


Co II 


1.160e-2 


M03 


2515.0725 


Si I 


2.110e-l 


M03 


1560.31 


CI 


7.740e-2 


M03 


2522.1238 


Col 


2.290e-l 


M03 


1562.001 


Si I 


3.758e-l 


M01.M03 


2523.6083 


Fel 


2.030e-l 


M03 


1565.5255 


Fell 


0.766e-5 


M03 


2568.7518 


All 


3.760e-2 


M03 


1573.92 


Co II 


3.700e-2 


M03 


2576.877 


Mn II 


3.610e-l 


M03 


1588.6876 


Fell 


1.480e-4 


M03 


2586.6500 


Fell 


6.910e-2 


M03 


1602.4863 


Ge II 


1.440e-l 


M99 


2594.499 


Mnll 


2.800e-l 


M03 


1608.4511 


Fell 


5.770e-2 


M03 


2600.1729 


Fell 


2.390e-l 


M03 


1610.194 


ScIII 


8.660e-2 


M03 


2606.462 


Mn II 


1.980e-l 


M03 


1611.2005 


Fell 


1.380e-3 


M03 


2653.2654 


All 


1.500e-2 


M03 


1631.170 


Si I 


2.051e-l 


M01.M03 


2669.948 


Al II] 


1.060e-5 


M03 


1651.164 


Mgl 


6.577e-4 


M01.M03 


2683.887 


VII 


1.090e-l 


M03 


1656.93 


CI 


1.490e-l 


M03 


2719.8329 


Fel 


1.220e-l 


M03 


1660.809 


O III]* 


1.590e-7 


M03 


2740.525 


VII 


1.010e-l 


M03 


1666.150 


O III]** 


2.090e-7 


M03 


2795.641 


Mnl 


5.650e-l 


M03 


1670.7886 


A1II 


1.740e00 


M03 


2799.094 


Mnl 


4.190e-l 


M03 


1674.5953 


PI 


3.250e-2 


M03 


2801.907 


Mnl 


2.900e-l 


M03 


1679.6969 


PI 


5.360e-2 


M03 


2796.3543 


Mg II 


6.155e-l 


M03 


1683.4116 


Mgl 


2.490e-3 


M03 


2803.5315 


Mg II 


3.058e-l 


M03 


1703.4119 


Nill 


6.000e-3 


M03 


2852.9631 


Mgl 


1.830e00 


M03 


1709.6042 


Nill 


3.240e-2 


M03 


2937.7623 


Fel 


1.810e-2 


M03 


1741.5531 


Nill 


4.270e-2 


M03 


2967.7646 


Fel 


4.380e-2 


M03 


1747.7937 


Mgl 


9.080e-3 


M03 


2984.4402 


Fel 


2.900e-2 


M03 


1751.9157 


Nill 


2.770e-2 


M03 


3021.5187 


Fel 


1.040e-l 


M03 


1753.8406 


Mgl 




M03 


3067.2379 


Till 


4.890e-2 


M03 


1774.9487 


PI 


1.690e-l 


M03 


3073.8633 


Till 


1.210e-l 


M03 


1782.8291 


PI 


1.130e-l 


M03 


3083.0462 


All 


1.770e-l 


M03 


1787.6481 


PI 


6.040e-2 


M03 


3131.54 


Bell 


4.981e-l 


M03 


1804.473 


Nill 


3.800e-3 


W99 


3187.372 


Til 


1.950e-l 


M03 


1807.3113 


S I 


9.050e-2 


M03 


3230.1221 


Ti II 


6.870e-2 


M03 


1808.0129 


Si II 


2.080e-3 


M03 


3242.9184 


Till 


2.320e-l 


M03 


1816.93 


Sill* 


1.790e-3 


M03 


3248.4739 


Cul 


4.340e-l 


M03 


1827.9351 


Mgl 


2.420e-2 


M03 


3274.8980 


Cul 


2.180e-l 


M03 


1845.5205 


Si I 


2.700e-l 


M03 


3303.529 


Nal 


1.345e-2 


M03 


1854.7184 


Aim 


5.590e-l 


M03 


3342.834 


Til 


1.740e-l 


M03 


1862.7910 


Aim 


2.780e-l 


M03 


3370.5336 


Nil 


2.430e-2 


M03 


1892.030 


Si III] 


2.690e-5 


M03 


3384.7304 


Till 


3.580e-l 


M03 


1901.773 


Fell 


7.000e-5 


M03 


3392.0165 


Nil 


0.999e-2 


M03 


1908.734 


CIII] 


1.687e-7 


M03 


3441.5918 


Fel 


2.360e-2 


M03 


1910.78 


Till 


2.020e-l 


M03 


3579.705 


CrI 


3.660e-l 


M03 


1941.2852 


Co II 


3.400e-2 


M03 


3581.947 


Sc II 


2.370e-l 


M03 


1977.5972 


Si I 


4.910e-2 


M03 


3594.507 


CrI 


2.910e-l 


M03 


2012.1664 


Co II 


3.680e-2 


M03 


3606.350 


CrI 


2.260e-l 


M03 
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Table A2(continued): Finding list for interstellar lines in spectra of QSOs: 1370-4500 A 



202(S IY70 


Zn II 


5 010e-l 


M03 


3636 498 

JUJU.T70 


Ti I 


2 520e-l 


M03 


7096 769 

Z.YJ . LAjy 


Cr II 


1 300e-3 


M03 


1641 877 


Sc II 


3.750e-l 


M03 


7026 4768 


Mo T 


1 130e-l 


M03 


1650 1476 


Fe I 


699e-4 


M03 


7056 756Q 


Cr II 


1 030e-l 


M03 


1680 Q61 1 


Fe I 


7 800p-1 


M03 


7067 7161 

_WU_ ._ jll 1 


Cr II 


7 590e-2 


M03 


1770 QQ78 

J / Lxi.yy Z.O 


Fe I 


4 11 Oe-2 


M03 


7067 6604 


Zn II 


7 460p-1 


M03 


1875 5788 


Fe I 


4 810p-1 


M03 


7066 1640 


Cr II 


5 170e-2 

J . 1 Z.VJC z. 


M03 


1861 00*58 


Fe I 


Z,. 1 / wc z. 


M03 


71 18 CT7 
^1 JO. oz / 


V II 


e 7 00p- 7 


M03 


1875 705 

DO 1 J. 1 UJ 


rw RSXSO 

DJAJU 


3 3R0e-7 


BvD 


7 HQ 7477 


Zn I 


1 47 Op 00 


M03 


187Q 871 
Jo / y .0 / J 




1.070e-3 


GvDB 


71 1Q 681 

Z, 1 J7.UOJ 


N II]* 


<S 860p-8 


M03 


1887 510 


VII DJATVI 


3 900e-3 


GvDB 


2143.450 


Nil]** 


8.190e-8 


M03 


3891.319 


CH BSXP0 


2.130e-3 


GvDB 


2167.4534 


Fel 


1.500e-l 


M03 


3934.7750 


Ca II 


6.267e-l 


M03 


2208.6666 


Si I 


5.750e-2 


M03 


3945.1224 


All 


1.170e-l 


M03 


2234.4472 


Fe II 


2.520e-5 


M03 


3958.812 


CH+APXS1 


3.400e-3 


GvDB 


2249.8768 


Fe II 


1.820e-3 


M03 


3969.5901 


Call 


3.116e-l 


M03 


2260.7805 


Fe II 


2.440e-3 


M03 


3982.888 


Til 


1.050e-l 


M03 


2264.1647 


All 


8.920e-2 


M03 


4045.2847 


KI 


5.680e-3 


M03 


2290.6934 


Nil 


1.280e-l 


M03 


4048.3565 


KI 


2.640e-3 


M03 


2299.663 


Mnll 


4.810e-4 


M03 


4227.918 


Cal 


1.770e00 


M03 


2305.714 


Mnll 


1.150e-3 


M03 


4233.740 


CH+APXS0 


5.500e-3 


GvDB 


2311.6723 


Nil 


3.765e-l 


M01.M03 


4301.5229 


CH ADXP0 


5.060e-3 


GvDB 


2320.7468 


Nil 


6.850e-l 


M03 


4430.12 


DB 4428 


l.OOOeOO 


dib 



References: BvD: Black, J. H. & van Dishoeck, E. F. 1988, ApJ, 331, 986 
GvBD: Gredel, R., van Dishoeck, E. F, & Black, J. H. 1993, A&A, 269, 477 
M91: Morton, D.C. 1991, ApJS, 77, 119; M03: Morton, D.C. 2003, ApJS, 149, 205 
MN: Morton, D. C. & Noreau 1994, ApJS, 95, 301 
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Tabic A3: Equivalent widths of weak lines for derivation of abundances 



Sub-sample Number 


3 


4 


24 


9 


n a 


13 a 


12 a 


10 a 


26 


Number of systems 


135 


133 


698 


404 


398 


405 


411 


405 


97 




Criteria 




w Mg II 
0.9-1.2 


W Mg II 
1.2-1.5 


a b 
A 


z abs 


1 





% 


z abs 


W Mg II 












< 0.2 


<1.3 


<19.1 


< 0.1 


> 19.1 


> 1.3 


>2.5 
A < 0.2 




E(B - V) 




<0.001 


< 0.001 


0.002 


0.006 


0.009 


0.011 


0.011 


0.012 


0.012 




Species 


/ 










Equivalent Width 


in mA 








1656.9 


C I 


0.1490 


<48 


<78 


16±4 


NA 


<24 


<40 


<39 


39±13 


49±18 


3303.6 


Na I 


0.0130 


<29 


<24 


< 13 


< 15 


< 15 


< 19 


<32 


<33 


<33 


2853.0 


Mg I 


1.8300 


171 ±22 


281 ± 16 


258±4 


341 ±20 


298 ± 13 


279±15 


379±21 


3O0±17 


596±16 


2026.5 


Mgl 


0.1130 


<23<>7) 


<39(> 1 1) 


< 14(> 10) 


<7(>14) 


<21(> 12) 


<29(> 1 1) 


<37(> 15) 


< 15(> 12) 


<82(>24) 


1827.9 


Mgl 


n?77n 


< 17 


<21 


<8 


< 17 


<9 


< 10 


<21 


< 10 


<30 




'a! t 




<24 


<28 


<16 


< 13 


< 13 


<17 


<26 


<26 


<38 






1060 


<54 


<45 


<12 


<21 


<24 


<27 


<54 


<33 


<43 




1 


^700 


< 17 


<20 


<8 


< 17 


< 10 


< 10 


<20 


< 10 


<30 


2515.1 


Si 


0.2110 


< 18 


< 19 


<8 


< 1 1 


<9 


< 1 1 


<20 


< 12 


<29 


2.i 3. i.l 


Si III 


4.25E-6 


<19 


<17 


<8 


<8 


<8 


<10 


<19 


<1 1 


<26 


1 SIIS.O 


Si II 


0.0021 


63±18 


79±18 


62±5 


76 ±14 


59±7 


59±7 


83±24 


58±7 


187±19 


1807.3 


S I 


0.0905 


23±6 


<24 


29±3 


36 ±6 


30±3 


26 ±3 


29±8 


27 ±4 


95±10 


.>9 >I.K 


Ca II 


0.6270 


92±36 


< 147 


102±13 


98±26 


98 ±36 


199±84 


226±75 


NA 


253±56 


.. 969.6 


Ca II 


O.i 120 


<77 


< 106 


53±13 


102±26 


86 ±26 


156±51 


< 124 


NA 


< 67 


3643.8 


Sc 11 


0.3750 


<39 


<38 


< 16 


<21 


< 19 


<25 


<48 


<41 


<66 


3582.0 


Sc II 


0.2370 


<37 


<39 


<16 


< 17 


< 17 


<24 


<42 


<39 


<63 


3384.7 


Ti II 


0.3580 


<37 


<51 


26 ±5 


44±10 


34±12 


41 d= 12 


<44 


<39 


100±28 


3242.9 


Ti II 


0.2320 


<29 


<49 


< 16 


27±8 


29 ±7 


51±14 


<33 


<32 


78 = 19 




Ti II 


0.1980 


<24 


<21 


<9 


< 15 


<9 


< 10 


< 18 


< 10 


<26 


3073.9 


Ti II 


0.1210 


<24 


<28 


< 12 


17±4 


< 13 


< 18 


<25 


<24 


<38 




Till 


0.0690 


<25 


<31 


< 12 


< 14 


< 13 


<16 


<30 


<28 


<39 


268.1.9 


V II 


0.1100 


< 17 


< 18 


<8 


<8 


<8 


<9 


< 17 


< 13 


<29 


2740.5 


V II 


0. 1000 


<20 


<20 


<8 


<9 


<10 


< 1 1 


< 18 


< 14 


<27 


2056.3 




0.1030 


<23 


49 ±15 


35±5 


30±9 


23 ±6 


< 17 


45 ±15 


22±8 


93±16 


2062.2 


Cr II 


0.0760 


<23 


<36 


24 


<22 


< 17 


< 17 


<34 


< 18 


61 


2066.2 


Cr II 


0.0510 


<23 


<24 


12(4) 


< 14 


< 1 1 


< 17 


<23 


< 15 


28±9 


2576.9 


Mn II 


0.3610 


27±11 


61 ±14 


47 ±4 


61±8 


51±6 


43±7 


76±17 


53±9 


143±14 


2594.5 


Mn II 


0.2800 


39±14 


43±13 


42±4 


54±7 


54±9 


47 ±8 


60± 12 


42±9 


108±22 


2606.5 


Mn II 


0.1980 


<23 


24±11 


21 ±4 


24±7 


23±8 


15±7 


<23 


<14 


60±18 


2484.0 


Fel 


0.5440 


< 17 


<19 


<8 


< 11 


<9 


<11 


<21 


<13 


<28 


2523.6 


Fel 


0.2030 


<19 


<19 


<8 


<11 


<9 


<10 


<20 


<12 


<31 


2374.5 


Fell 


0.0313 


180±18 


281 ±15 


274 ±4 


307 ±7 


252 ±8 


264±9 


416±18 


253±11 


761 ±14 


2260.8 


Fell 


0.0024 


<22 


58±16 


43±4 


61±11 


48±7 


53±10 


59±17 


44±10 


137±16 


2249.9 


Fell 


0.0018 


<22 


36±11 


31 ±4 


42±8 


32±7 


30 ±8 


39±12 


19±9 


123±16 


2234.2 


Fell 


2.52E-5 


<18 


<18 


<8 


<10 


<8 


<10 


<17 


<12 


<27 


2367.6 


Fell 


2.16E-5 


<18 


20 ±6 


<9 


<11 


<9 


<10 


<20 


<12 


<26 


2012.2 


Co II 


0.0370 


< 18 


< 23 


13±4 


< 12 


< 9 


< 12 


< 19 


< 12 


< 26 


1941.3 


Co II 


0.0340 


<22 


<21 


<8 


<13 


<10 


<12 


<20 


<11 


<24 


1741.6 


Ni II 


0.0430 


<32 


<36 


25±4 


<31 


25±9 


29 ±7 


71 ±22 


30±9 


108±12 


1709.6 


Ni II 


0.0320 


<28 


<37 


20 ±6 


<25 


<14 


<14 


58±14 


<14 


62±13 


1751.9 


Ni II 


0.0280 


<22 


<32 


<10 


<26 


17±4 


23±5 


<23 


< 27 


40±9 


1703.4 


Ni II 


0.0060 


<30 


<37 


<12 


<25 


<15 


<15 


<24 


<14 


<39 


1804.5 


Ni II 


0.0038 


<54 


<54 


<14 


<42 


<21 


<21 


<72 


<21 


<58 


2139.3 


Znl 


1.4700 


<21 


<26 


<8 


<12 


<10 


<11 


<20 


<12 


<26 


2062.7 


Znll 


0.2460 


<23 


>35 


21 


>33 


>26 


<17 


>25 


>9 


30 


2026.3 


Znll 


bl c 


<23 


74 ±26 


35±4 


40±17 


21 ±6 


29 ±9 


62±19 


<15 


112±14 


2062.4 


Zn II 


bl c 


<23 


71 ±2 


45±4 


55±9 


43±10 


<17 


59±14 


28±8 


91±13 



a The numbers in the samples for measuring equivalent width were greater by 1,1, 1 and 2, for sub-samples 11, 13, 12 and 10, than what is listed here. The numbers shown for number of objects in 
the sub-sample and for the extinction of the sub-sample are from Table 1. 
k A is the observer frame colour excess A(g — i) 

c The equivalent widths for these lines are the actual measured values for the blended lines. 



Dust and abundances in QSO absorption line systems 



Table A3 (continued): Equivalent widths of weak lines for derivation of abundances 



Sample Number 




1 


16 


8 


17 


19 


25 


27 


Number of systems 




809 


369 


251 


85 


58 


111 


48 


Criteria 




W Mg II 


W Fe II^Mg II 


W Mg II 


W A 1 II™ M g I 


Fe II A 2260 


A >0.2 


W Mg II 






>0.3 


>0.58 


>2 


< 1.54 


present 




>2.5 
A>0.2 


E{B - V) 




0.013 


0.014 


0.0320 


0.034 


0.034 


0.081 


0.085 




X 


Speeies 






Equivalent Width in mA 








1656.9 


CI 


19±4 


46±13 


50± 12 


79±21 


< 57 


57±19 


114±32 


3303.6 


Nal 


<12 


<17 


<21 


<81 


<37 


31 


<67 


2853.0 


Mgl 


282±4 


397±10 


585±10 


691±35 


689±16 


536±15 


795 ±26 


2026.5 


Mgl 


<32(>11) 


<34(>16) 


<77(>23) 


<37(>28) 


<76(>28) 


>48(>21) 


<240(>32) 


1827.9 


Mgl 


<8 


<14 


<17 


<22 


<32 


<28 


<48 


3083.1 


All 


<14 


<14 


<24 


<46 


<35 


<34 


<58 


2367.8 


All 


< 11 


< 14 


<27 


<93 


<41 


<48 


<85 


1845.5 


Si I 


<8 


< 13 


< 17 


<22 


<31 


<29 


<49 


2515.1 


Si I 


<8 


< 12 


< 16 


<29 


<33 


<25 


<48 


2335.1 


Si III 


<7 


< 12 


<16 


<27 


<25 


<30 


<48 


1808.0 


Si II 


68±5 


121±10 


167±11 


145 ±22 


240±15 


175±24 


212±37 


1807.3 


SI 


32±3 


61 ±4 


87±5 


59±9 


124±11 


91±10 


115±18 


3934.8 


Ca II 


110±16 


254 ±64 


201 ±40 


NA 


254 ±47 


307 ±86 


752±139 


3969.6 


Ca II 


57±13 


156±42 


101 ±26 


NA 


141±36 


99±21 


< 163 


3643.8 


Se II 


<16 


<26 


<31 


< 106 


<50 


<59 


<82 


3582.0 


Se II 


< 16 


<24 


<30 


< 103 


<52 


<62 


<86 


3384.7 


Ti II 


31 ±5 


67±16 






160±16 


120±29 


193±56 


3242.9 


Ti II 


32±5 


49 ±10 


64±11 


<94 


120±16 


51 d= 13 


<432? 










<15 


<20 








3073.9 


Ti II 


<11 


<15 


<24 


<46 


<39 


<33 


<53 


3230.1 


Ti II 


<12 


<17 


<22 


<61 


<32 


<35 


< HI 


2683.9 


VII 


<8 


<11 


<17 


<29 


<29 


<28 


<42 


2740.5 


VII 


<8 


<14 


<18 


<29 


<27 


<27 


<43 


2056.3 


Cr II 


39±4 


72±11 


101 ±10 


89±17 


147±14 


77±13 


154±27 


2062.2 


Cr II 


40 


56 


78 


70 


123 


75 


132 


2066.2 


Cr II 


14±4 


39±10 


54±10 


50±16 


100±13 


72±18 


109 ±27 


2576.9 


Mn II 


54±4 


U4±ll 


140±9 


146±21 


231±14 


156±17 


237 ±26 


2594.5 


Mn II 


50±4 


78±10 


129±11 


103±16 


167±12 


162±19 


277 ±32 


2606.5 


Mn II 


26 ±4 


49 ±9 


79 ±10 


64±14 


142±13 


96±18 


147 ±30 


2484.0 


Fel 


<8 


<13 


<18 


<26 


<33 


<28 


<50 


2523.6 


Fel 


<8 


<12 


<16 


<30 


<34 


<25 


<48 


2374.5 


Fell 


298±4 


503 ±12 


702 ±9 


476±31 


894±14 


624±16 


958±28 


2260.8 


Fell 


46 ±4 


99±11 


113±9 


125 ±23 


280±10 


90±14 


145 ±18 


2249.9 


Fell 


35±4 


63 ±10 


105±10 


58±17 


163±15 


113±20 


215±32 


2234.2 


Fell 


<8 


<12 


<17 


<27 


<24 


<30 


<51 


2367.6 


Fell 


<8 


<13 


<16 


<27 


<25 


< 30 


<50 


2012.2 


Co II 


15±4 


< 11 


56±13 


< 20 


57±17 


< 30 


< 47 


1941.3 


Co II 


<8 


<10 


<16 


<19 


<27 


<26 


<43 


1741.6 


Ni II 


28±4 


84±24 


90±8 


80 ±20 


124±16 


87±19 


100±23 


1709.6 


Ni II 


16 ±4 


22±5 


83±17 


<29 


117±24 


<41 


<63 


1751.9 


Nin 


<10 


<16 


51 ±17 


52±13 


<40 


76 ±20 


93 ±30 


1703.4 


Ni II 


<11 


<18 


<38 


<28 


<45 


<42 


<66 


1804.5 


Ni II 


<14 


<30 


<32 


<66 


<36 


<72 


<110 


2139.3 


Znl 


<8 


<12 


<17 


<25 


<28 


<25 


<47 


2062.7 


Znll 


8 


21 


36 


31 


102 


47 


26 


2026.2 


Znll 


40±4 


55±8 


111±8 


68±127 


176±14 


95±15 


266 ±29 


2062.7 


Znll 


48±4 


77±10 


114±8 


101 d= 15 


225 ±14 


122±14 


158±20 
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Table A4: Column densities of weak lines 



Sample Number 


3 


4 


24 


9 


11 


13 


12 


10 


Number of systems 


135 


133 


698 


404 


398 


405 


411 


405 


Criteria 


W Mg II 


W Mg II 


A a 


z abs 


i 





i 


z abs 




0.9-1.2 


1.2-1.5 


<0.2 


<1.3 


< 19.1 


<0.1 


> 19.1 


> 1.3 


E(B — V) 


<0.001 


<0.001 


0.002 


0.006 


0.009 


0.011 


0.011 


0.012 


Species 








Column densities' 5 in cm 1 








N c /10 20 
H I 


<0.5 


<0.5 


0.9 


2.8 


4.2 


5.0 


5.0 


5.4 


log N H i 


<19.7 


<19.7 


20.0 


20.4 


20.6 


20.7 


20.7 


20.7 


Mgll 


14.2a[0.5] 


14.9[1] 


14.7a[l] 


14.9[1] 


14.9[1] 


14.9[1] 


15.0[0.9] 


14.8[0.9] 


Al II 


13.1a[0.2] 


13.6a[0.5] 


13.5a[0.5] 


13.6b[0.5] 


13.5a[0.5] 


13.5a[0.5] 


13.6a[0.5] 


13.6a[0.5] 


Aim 


13.0b 


13.1a 


13.1a 


13.2a 


13.1a 


13.1a 


13.3a 


13.1a 


Si II A 1808 


15.1c 


15.1b 


15.0a 


15.1b 


15.0a 


15.0a 


15.0a 


15.0a 


Si II A1526 d 


14.0a 


14.1b 


14.1a 


NA 


14.1e 


14.1a 


14.2a 


14.1a 


Ca II 


12.0c 


< 12.0 


<122 


12.4b 


12.3b 


12.6b 


12.4c 


NA 


Sc II 


<12 


<11.9 


< 1 1.6 


<11.7 


<11.6 


< 1 1.8 


< 12.0 


< 12.0 


Ti II 


< 1 1.9 


<11.9 


11.9b 


12.1b 


12.1b 


12.4c 


<12.0 


< 1 1.9 


VII 


<12.4 


<12.4 


<12.1 


<12.1 


<12.1 


<12.1 


< 12.4 


< 12.3 


CrII 


<12.7 


13.1c 


12.9b 


12.9b 


12.8b 


<12.6 


13.1c 


12.8c 


Mn II 


12.1c 


12.5c 


12.3a 


12.5a 


12.4a 


12.4b 


12.6b 


12.4b 


Fell 


<14.2 


14.7b 


14.7b 


14.8b 


14.7b 


14.7b 


14.7b 


14.6b 


Coll 


< 13.3 


<13.3 


< 12.9 


<13.1 


< 13.0 


<13.1 


< 13.3 


< 13 


Ni II 


<13.4 


<13.4 


(13.4) 


< 13.4 


13.3c 


13.4b 


13.8c 


13.4c 


Znll 


< 12.3 


> 12.6 


12.4 


> 12.4 


>12.4 


<12.3 


>12.4 


>12.0 


CI 


<13 


<13 


12.6c 


NA 


< 12.8 


<12.7 


<13 


13.0c 


Mgl 


12.1-12.6 


12.4-12.8 


12.3-12.5 


12.4-12.6 


12.4-12.6 


12.4-12.7 


12.5-12.8 


12.4b 


All 


< 12.2 


<122 


<11.9 


<11.9 


<12.0 


<12.1 


<122 


<12.3 


Si I 


< 12.2 


<122 


<11.8 


<12.0 


<12.0 


<12.0 


<122 


<12.0 


Fel 


< 1 1.8 


<11.8 


<1 15 


<11.6 


<11.5 


<11.6 


<11.9 


< 1 1.6 



a A is the observer frame colour excess A(g — i) 

b Letters (lower case) indicate errors, dex, a=0-0.049, b=0.05-0.099, c=0.1-0.149, etc. () indicates uncertainties or inconsistent raw data. 

C i comes from using the SMC extinction curve to derive E(B — V) and the using the SMC ratio of N H j to E(B — V) from Welty et al. 2005 (which is 
similar to previous values). 

d We used the column densities from Si II A 1808. The values from Si II A 1526 are to show that the saturation corrections for Si II are consistent with the saturation 
corrections for Al II and Mg II derived from Fe II (in brackets for the entries for these two species.) 



Table A4 (Continued): Column densities of weak lines 



Sample Number 


26 


1 


16 


8 


17 


19 


25 


27 


Number of systems 


97 


809 


369 


251 


85 


58 




48 


Criteria 


w Mg II 


WMg II 


WFo Il'WMg II 


WMg II 


W A 1 IlAV M g i 


FeIIA2260 


A > 0.2 


w Mg II 




>2.5 


>0.3 


>0.58 


>2.0 


<1.54 


present 




>2.5 




A <02 














A >02 


E(B - V) 


0.012 


0.013 


0.014 


0.032 


0.034 


0.034 


0.081 


0.085 


Species 








Column densities in cm 1 








N H l/10 20 


5.4 


5.9 


6.3 


14.4 


15.3 


15.3 


36.0 


38.0 


log N H I 


20.7 


20.8 


20.8 


21.2 


21.2 


21.2 


21.6 


21.6 


Mg II 


15.111.0] 


14.7[0.9] 


14.9[0.9] 


15.0a[0.9] 


15.1[1.0] 


15.2[1.2] 


14.9[0.9] 


15.2[1.0] 


Al II 


14.0a[0.6] 


13.5a[0.5] 


13.7a[0.5] 


13.9a[0.5] 


13.9a[0.5] 


14.1a[0.7] 


13.8[0.5] 


14.0a[0.5] 


A1III 


13.4a 


13.1a 


13.3a 


13.4a 


13.3a 


13.3a 


13.5a 


13.6a 


SiIIA1808 


15.5a 


15.1a 


15.3a 


15.4a 


15.4b 


15.6a 


15.5a 


15.6b 


SiIIA1526 


14.5a 


14.1a 


14.4a 


14.5a 


14.4a 


14.5a 


14.6a 


14.8a 


Call 


>12.5 


12.1 


12.6c 


12.4c 


NA 


12.6b 


12.6c 


>12.6 


Sc II 


< 12.2 


< 1 1.6 


<11.8 


<11.9 


< 12.4 


< 12.1 


<12.1c 


<12.3 


Ti II 


12.5c 


11.9b 


12.4b 


12.4b 


12.6e 


12.7b 


12.3c 


< 12.5 


VII 


<12.6 


<12.1 


< 12.2 


< 12.4 


< 12.6 


< 12.6 


<12.6 


<12.8 


CrII 


13.3b 


13.0b 


13.3b 


13.3b 


13.4c 


13.7b 


> 13.6c 


13.8b 


Mn II 


12.9b 


12.4a 


12.7a 


12.9a 


12.8a 


> 13.1 


13.0a 


> 13.2 


Fell 


15.1b 


14.7a 


15.0a 


15.1b 


15.1b 


15.4a 


>14.9 


> 15.3c 


Coll 


< 13.3 


< 12.8 


<13.0 


< 13.1 


<13.2 


<13.4 


< 13.4 


< 13.6 


Ni II 


(13.9c) 


(13.4c) 


13.9c 


<(13.9c) 


13.9b 


(13.9c) 


> 14.1 


(>14.2) 


Znll 


12.5 


11.9 


12.4 


12.6 


12.5 


13 


> 12.7 


> 12.4 


CI 


13. le 


12.7c 


13.1b 


13.1c 


13.4b 


< 13.2 


13.6: 


13.5d 


Mgl 


12.7-13.1 


12.3-12.7 


12.5-12.8 


12.6-13.1 


12.5-12.7 


12.7-13.1 


12.6-12.9 


12.9-13.8 


All 


< 12.4 


< 1 1.9 


< 12.0 


< 12.2 


< 12.5 


< 12.4 


< 12.4 


<12.6 


Si I 


< 12.4 


<11.8 


<12.0 


< 12.1 


< 12.4 


< 12.5 


< 12.3 


< 12.6 


Fel 


<12.0 


< 1 1.4 


<11.7 


<11.8 


<12.0 


< 12.1 


<12.0 


<12.2 



